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We report a measurement of the differential cross section of 7 r° pair production in single-tag 
two-photon collisions, 7*7 —> 7 r° 7 r°, in e + e _ scattering. The cross section is measured for Q 2 
up to 30 GeV 2 , where Q 2 is the negative of the invariant mass squared of the tagged photon, 
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in the kinematic range 0.5 GeV < W < 2.1 GeV and | cos (9* | < 1.0 for the total energy and 
pion scattering angle, respectively, in the 7*7 center-of-mass system. The results are based on a 
data sample of 759 ftU 1 collected with the Belle detector at the KEKB asymmetric-energy e + e“ 
collider. The transition form factor of the /o(980) and that of the /2(1270) with the helicity-0, -1, 
and -2 components separately are measured for the first time and are compared with theoretical 
calculations. 

PACS numbers: 12.38.Qk, 13.40.Gp, 14.40.Be, 13.60.Le, 13.66. Be 


I. INTRODUCTION 

Single-tag two-photon production of a C-even 
hadronic system (A4), 7*7 —>• A4, is an important re¬ 
action to investigate the nature of strong interactions 
in the low energy region, where perturbative Quan¬ 
tum Chromodynamics (QCD) cannot be applied. It 
also provides valuable information on the Q 2 depen¬ 
dence of the transition form factor (TFF), where Q 2 
is the negative of the invariant mass squared of the 
tagged photon. This reaction can be studied through 
the process of e + e~ —> e ± (e =F )A / l, where (e^) in¬ 
dicates an undetected electron or positron, and re¬ 
sults of the measurement can be directly compared 
to QCD-based theoretical predictions. Diehl, Gous- 
set and Pire considered this process at large Q 2 and 
small W (< 1 GeV) in terms of constituent-hard 
scattering and generalized distribution amplitudes 
and predicted a sizable cross section at LEP and 
B factories [I]; it is indeed the case at a B factory 
as reported here. Based on this framework, Braun 
and Kivel pointed out that the measurement of the 
TFF of the /2(1270) will be useful to cleanly deter¬ 
mine a gluon admixture in tensor mesons at large 
enough Q 2 [2j. In addition, a data-driven disper¬ 
sive approach was suggested recently [ 3 ], allowing a 
more precise estimate of the hadronic light-by-light 
contribution to the anomalous magnetic moment of 
the muon (g — 2). 

Schuler, Berends, and van Gulik, who had calcu¬ 
lated meson TFFs based on a heavy quark approxi¬ 
mation [1] , found that their calculations apply well to 
light mesons as well with only minor modifications. 
The predicted Q 2 dependence of TFFs for mesons 
with J PC = 0 ++ and 2 ++ is summarized in Table [I] 
where the 7*7 center-of-mass (c.m.) energy W is re¬ 
placed by the resonance mass M and A represents 
the total helicity of the two incident photons. Note 
that the helicity-1 (A = 1) state is allowed when a 
photon is off the mass shell. According to Table [IJ 
TFFs for the helicity-0 and -1 components of a tensor 
meson grow with Q 2 , a prediction which is amenable 
to investigation. 


TABLE I: Predicted Q 2 dependence of transition form 
factors of mesons for various helicities A of two incident 
photons [4]. Each term has a common divisor of (1 + 
Q 2 /M 2 ) 2 . 


J PC 

Q 2 dependence ^-r 

(i + Af) 


A = 0 

A = 1 

A = 2 

0 ++ 

( 1 + Sa? j 

- 

- 

2 ++ 

Q 2 

v'Sm 2 

Vo 2 

V2M 

1 




Recently, Pascalutsa, Pauk, and Vanderhaeghen 
have formulated sum rules for 7*7* fusion cross sec¬ 
tions, finding several new exact superconvergence re¬ 
lations that are integrated to zero |S]. They derive 
two predictions for the helicity-2 TFF of the /2(1270) 
from two sum rules in the case of one virtual and 
one real photon under the assumption that the sum 
rules are saturated by low mass resonances includ¬ 
ing the /2(1270). In one sum rule, the integrand has 
contributions from pseudoscalar mesons and tensor 
mesons. In the other, axial-vector mesons and ten¬ 
sor mesons contribute to its integrand. The first 
(second) sum rule gives the helicity-2 TFF of the 
/2(1270) in terms of TFF information of the rj and 
7 / (/i(1285) and /i(1420)). With pseudoscalar (P) 
mesons, the helicity-2 TFF of the /9(1270) is given 

by 



where A,, and are the pole masses and / = 
c v/( c v + <v)> with cp = r 77 (P)/m.p. The relevant 
parameters are summarized in Table [TlJ In another 
sum rule for axial-vector (A) mesons, the helicity-2 
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TABLE II: Parameters of the 7 , rf [Q, /i(1285), and 
/i(1420) [3 [§]. r 77 for an axial-vector meson shall read 
f 77 defined in Eq. |3j|. 


Meson 

M M (MeV/c 2 ) 

r 77 (keV) 

A M (MeV/c 2 ) 

V 

547.853 ±0.024 0.510 ±0.026 

774 ± 29 

v 1 

957.78 ± 0.06 

4.29 ±0.14 

859 ± 28 

/i(1285) 

1281.8 ± 0.6 

3.5 ±0.8 

1040 ± 78 

/i(1420) 

1426.4 ±0.9 

3.2 ±0.9 

926 ± 78 



TFF of the /2(1270) is given by 


FfAQ 2 ) 



where A f 1 and are the pole masses and / = 
Cfi/( c f 1 +cy'), with ca = 3r 77 (A)/m^. The effec¬ 
tive two-photon width of the A resonance is defined 
as 

M 2 

F 77 (A) = Jim o g#r(A -± 727 t), (3) 


FIG. 1: Feynman diagram for the process e + e“ —> 
e(e)n°7r°. p i, p), p 2 , and p' 2 are the four-momenta of 
the incident and scattered electron or positron, q\ , q 2 are 
those of the virtual and real photons and ki, k 2 are those 
of the produced 7r° mesons. 




where T(A —► 7 |, 7 r) is the parameter of the axial- 
vector meson A decaying into a virtual longitudinal 
photon and a real transverse photon. The relevant 
parameters are also summarized in Table [TT] 


FIG. 2: Definition of the coordinate system for 7*7 —> 
7r 0 7T°. The incident 7 * and 7 are along the z* axis, the 
tagged e ± is in the x*z* plane with p x *tag > 0, and a 7 r° 
is produced at angles (6*, ip*). 


Experimentally, for pseudoscalar mesons, the TFF 
of the 7 r° meson has been measured recently by 
BaBar j9] and by Belle [10], and those of 7 and r/m 
and rj c | 12 | by BaBar for Q 2 < 40 GeV“. 

Two-photon production of axial-vector mesons, 
/i(1285) and /i(1420), which is interpreted as a 
two-photon fusion of a longitudinal (helicity- 0 ) and 
a real photon, was studied by the L3 collabora¬ 
tion, who measured the parameters listed in Ta¬ 
ble |IT] |7J E]. For scalar or tensor mesons, no signifi¬ 
cant data for the high-Q 2 region beyond the p-meson 
mass scale exist to be compared with QCD predic¬ 
tions; only yields consistent with zero were reported 
for 7*7 —► / 2 (1270) —► 7 r + 7 T - at Q 2 > 1.0 GeV 2 by 
the TPC/Two-Gamma collaboration m- 

We report a measurement of the process e + e — ► 
e(e)7r°7r°, where one of the e ± is detected together 
with 7r°7T 0 while the other is scattered in the for¬ 
ward direction and undetected. The Feynman dia¬ 
gram for the process is shown in Fig. |T] where the 


four-momenta of the particles are defined. We con¬ 
sider the process 7*7 —> Tr°n° in the c.m. system 
of the 7 * 7 . We define the x*j/*z*-coordinate system 
as shown in Fig. [ 2 ] at fixed values of W and Q 2 ; the 
asterisks here denote the coordinate system that is 
used for angular variables for the differential cross 
sections. One of the 7 r° mesons is scattered at an¬ 
gles (0*, ip*). Because of the identical particles in 
the final state and P symmetry in the reaction, only 
the region where 9* < 7r/2 and 0 < ip* < it is mean¬ 
ingful. The 2 *-axis is along the incident 7 * and the 
x* z* plane is defined by the tagged such that 
Px* tag > 0, where p tag is the 3-momentum of the 
tagged e ± . 

The differential cross section for 7*7 —> n 0 n 0 is 
given by [14] 


d<i( 7*7 —> n°Tr°) 
d£l 


2 

Y. tn cos(n<p*), 

n=0 


( 4 ) 
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with 


to — \M ++ \ 2 + |M_|—1~ + 2eo|Mo+|“, (5) 

h = 2 ei a?(( m;_-m; + )m 0+ ), (6) 

i 2 = -2e 0 5R(M;_M ++ ). (7) 


Here, M ++ , Mq + , and M_|_are helicity amplitudes 

whose subscripts +, —, and 0 indicate the helicity 
state of the incident virtual photon along, opposite, 
and transverse to the quantization axis, respectively, 
and eg and ei are given by 


eo 


ei 


1 — x 

1 ! T 5 " * 

1-x + f- 


(2-x)y/^ 

1-^+f 


( 8 ) 

( 9 ) 


In turn, x is defined as 


(gi • 92) 

(Pi - 92 )’ 


( 10 ) 


where gi, g 2 ,pi, and p 2 are the four-momenta of the 
virtual and real photons and the incident electron 
and positron, respectively, as defined in Fig.[1] When 
Eq. Q is integrated over ip*, we obtain 


da{ 7*7 —> 7r°7T 0 ) 

47tg?| cos 9* | 


— \M ++ \ 2 + |M_|—1“ + 2eo|Mo+|‘ 


( 11 ) 

The total cross section is obtained by integrating 
Eq. (Ill over cos0*. It can be written as 


Otot(7*7 7r ' ^) = a TT + eoVLT, (12) 


where cttt {ctlt) corresponds to the total cross sec¬ 
tion of two photons, both of which are transversely 
polarized (one transversely and the other longitudi¬ 
nally polarized); as Q 2 —> 0, the second term van¬ 
ishes and ott approaches the total cross section of 
real photon-photon scattering. 

Neutral-pion pair production in the final state 
e(e)7r°7T 0 is different from the corresponding charged- 
pair process, e(e)7r + 7r~: the 7r°7r° is a pure (7-even 
state, whereas the 7 r + 7 r _ is a mixture of (7-even and 
(7-odd states. Thus, the 7 r° 7 r° state has no contribu¬ 
tion from single-photon production (“bremsstrahlung 
process”), whose effect must be considered in two- 
photon production of 7r + 7r - . 

In this paper, we report for the first time a mea¬ 
surement of the cross section for the process 7*7 —► 
7r°7T 0 up to Q 2 = 30 GeV 2 , from which we extract the 
TFF of the / o (980) and helicity-0, -1, and -2 TFFs 
of the / 2 (1270). 


This article is organized as follows. Section |TT] 
briefly describes the Belle detector and the data 
sample used in this measurement. The Monte 
Carlo (MC) program used for producing simulated 
events and for efficiency determination is described 
in Sec. |III| Selection of events and comparison with 
MC data are explained in Sec. |IV| Section [V] is de¬ 
voted to estimation of possible backgrounds. The 
differential cross section is derived and its systematic 
uncertainties are estimated in Sec. [VI] In Sec. |VII| 
the cross section is parameterized and fitted to ex¬ 
tract the TFFs of the /o(980) and the helicity-0, -1, 
and -2 components of the / 2 ( 1270) as a function of 
Q 2 , which are compared to theoretical predictions. 
Finally, Sec. |VIlT| provides the summary and conclu¬ 
sion. 


II. EXPERIMENTAL APPARATUS AND 
DATA SAMPLE 

In this section, we briefly describe the Belle de¬ 
tector and the data sample. We use a 759 fb -1 
data sample recorded with the Belle detector dans 
at the KERB asymmetric-energy e + e _ collider m- 
We combine data samples collected at several beam 
energies: at the Y(4 S) resonance {yfs = 10.58 GeV), 
where the beam energy for the electron (positron) 
beam is 8 GeV (3.5 GeV), and 60 MeV below it 
(637 fb -1 in total); at the T(35) resonance (y/s = 
10.36 GeV, 3.2 fb -1 ); and near the T(5S) resonance 
(y/s = 10.88 GeV, 119 fb -1 ). Correspondingly, when 
combining the data, the slight dependence of the 
two-photon cross section on beam energy is taken 
into account as described in Sec. ED 

This analysis is performed in the “single-tag” 
mode, where either the recoil electron or positron 
(hereinafter referred to as an electron) alone is de¬ 
tected. As described in Sec. |IV| in more detail, 
we restrict the virtuality ( Q 2 ) of the untagged-side 
photon to be small by imposing a strict transverse- 
momentum balance between the tagged electron and 
the final-state neutral pion pair in the e“*"e - c.m. 
frame with respect to the beam axis. In this article, 
we refer to events tagged by an e + or an e~ as “p-tag” 
(positron-tag) or “e-tag” (electron-tag), respectively. 


A. Belle detector 

A comprehensive description of the Belle detec¬ 
tor is given elsewhere [Dp fl 6 ] . In the following, we 
describe only the detector components essential for 
this measurement. Charged tracks are reconstructed 
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from the drift-time information in a central drift 
chamber (CDC) located in a uniform 1.5 T solenoidal 
magnetic field. The z axis of the detector and the 
solenoid is along the positron beam direction, with 
the positron beam pointing in the —z direction. The 
electron-beam direction at the collision point is 22 
mrad from the z axis. The CDC measures the longi¬ 
tudinal and transverse momentum components, i.e., 
along the 2 axis and in the rip plane perpendicular 
to the beam, respectively. Track trajectory coordi¬ 
nates near the collision point are provided by a sil¬ 
icon vertex detector. Photon detection and energy 
measurements are performed with a CsI(Tl) electro¬ 
magnetic calorimeter (ECL) by clustering of the en¬ 
ergy deposits in the crystals of the electromagnetic 
shower from an energetic photon, where a possible 
connection of a charged track to the cluster is exam¬ 
ined by extrapolating the track trajectory. Electron 
identification (ID) is based on E/p, the ratio of the 
calorimeter energy to the track momentum. 


B. Data sample 


To be recorded, events of interest here must satisfy 
one of the two ECL-based triggers: the HiE (High- 
energy threshold) trigger and the Clst4 (four-energy- 
cluster) trigger [T5] . 

The HiE trigger requires that the sum of the en¬ 
ergies measured by the ECL in an event exceed 
1.15 GeV but that the event be not Bhabha-like; 
the latter requirement is enforced by the absence 
of the CsiBB trigger (“Bhabha-veto”), which is de¬ 
signed to identify back-to-back Bhabha events m- 
For the purpose of monitoring trigger performance, 
we record one in 50 events that satisfy the CsiBB 
trigger (i.e., prescaled by a factor of 50). 

The Clst4 trigger requires at least four energy clus¬ 
ters in the ECL with each cluster energy larger than 
0.11 GeV. This trigger is not vetoed by the CsiBB 
because the Bhabha-event rate is manageable for the 
Clst4 trigger sample. Five clusters are expected in 
total in the signal events of interest if all the final- 
state particles are detected in the region where the 
ECL trigger is sensitive. 

We do not use information from the charged-track 
triggers because they require two or more charged 
tracks whereas our signal has only one. 


III. SIGNAL MONTE CARLO CODE 
A. Signal Monte Carlo, TREPSBSS 


We use the signal Monte Carlo (MC) generator 
TREPSBSS, which has been developed to calcu¬ 
late the efficiency for single-tag two-photon events, 
e + e _ -A e(e)X, as well as the two-photon luminosity 
function for 7*7 collisions at an e + e _ collider. 

TREPSBSS implements Eqs. (2.16) to (2.20) of 
Ref. |19| and is based on the MC code in Ref. | 20 | . 
which was modified to match the single-tag config¬ 
uration. We regard Eq. (12 1 as the total cross sec¬ 
tion of the 7*7 collisions, according to Eq. (2.16) 
of Ref. (IS], although eo is a variable depending on 
experimental conditions. It is possible to estimate 
e 0 for specific experimental conditions by taking the 
average value of eo calculated for the selected sig¬ 
nal events in MC or experimental data. Under our 
experimental conditions, eo ranges from 0.7 to 0.9. 
However, ctlt cannot be separated from er TT based 
on this information only. 

In TREPSBSS, the following kinematical variables 
are used for characterizing 7*7 collisions of the gen¬ 
erated events and for an integration to calculate the 
two-photon luminosity function: Q\ and Q\ (the ab¬ 
solute value of momentum transfer squared of the 
highly virtual and the less virtual incident photons, 
respectively), W (c.m. energy of the incident 7*7 
system), w 2 (the energy of the photon with the 
smaller virtuality), and A ip (the azimuthal-angle dif¬ 
ference between the two virtual photons). The choice 
of kinematical variables is discussed in Ref. [21]. 

We always retain the condition Q\ > Q\ for the 
virtuality of the two colliding photons by requiring 
Qi > 3.0 GeV 2 and Q\ < 1.0 GeV 2 for the ranges of 
integration, event generation, and selection. This Q\ 
range is sufficient to generate signal events over the 
kinematical region for the \Yip* t | selection criterion. 

The form-factor effect for the photon with the 
smaller virtuality is assumed to follow the 1/(1 + 
Q|/to 2 ) 2 dependence, where m p is the p-meson 
mass, 0.77 GeV/c 2 . The cross section defined be¬ 
low is that extrapolated to Q\ = 0, assuming this 
Q 2 dependence. 

The distribution of Q\ (= Q 2 ) in the MC sample is 
arbitrary and should not affect the final result. Con¬ 
ventionally, we choose the distribution corresponding 
to the flat form factor in order to retain high statis¬ 
tics of the signal MC in the high-Q 2 region but we 
weight the MC sample by an additional 1 /Q 2 factor 
to model a more realistic distribution for the effi¬ 
ciency derivation in each Q 2 region discussed in this 
analysis, and for comparison of distributions between 
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the signal-MC samples and the experimental data. 

The luminosity function, which is a conversion fac¬ 
tor between the e + e _ -based cross section a ee and 
the 7*7-based cross section <t 7 » 7 , is calculated in the 
same code. It is defined by 


d 2 c 


rj-T * 
_ 2 u -^7 7 


dQ 2 dW dQ 2 dW 


oy .7 (Q 2 , WO > (!3) 


where 2L 7 » 7 corresponds to the value of e + e _ -based 
integrated cross section per unit 7*7-based cross sec¬ 
tion in Eq. (2.16) of Ref. [15], and the differential lu¬ 
minosity function d 2 L 1 - 1 /dQ 2 dW is calculated by 
performing a three-dimensional numerical integra¬ 
tion over Q 2 , Cl> 2, and Aip. The factor of two incor¬ 
porates for the contributions from the e-tag and p- 
tag processes included in the entire e + e _ -based cross 
section. 


momentum of the ISR photon to determine the true 
Q 2 value in signal-MC events and to study the Q 2 de¬ 
pendence of the detection efficiency. We correct the 
experimental Q 2 dependence for the ISR effect using 
factors obtained from the signal MC (see Sec. VID1. 
We choose 15 different W points between 0.4 GeV 
and 2.5 GeV for the calculation of the luminosity 
function and event generation. 

We use a GEANT3-based detector simulation [51] 
to study the propagation of the generated particles 
through the detector; the same code as for the exper¬ 
imental data is used for reconstruction and selection 
of the MC simulated events. We thus obtain the se¬ 
lection efficiency as functions of Q 2 , W, and | cos0*| 
for a flat ip* dependence of the differential cross sec¬ 
tion. A correction for the observed <p* dependence is 
discussed in Sec. I VI Cl 


IV. EVENT SELECTION 

B. Event generation 


The event generation is performed using the same 
integrand but including initial-state radiation (ISR) 
effects from the tag-side electron. Inclusion of ISR 
changes the kinematics and Q 2 significantly. Mean¬ 
while, ISR from the untagged side has little effect 
because an ISR photon is nearly parallel not only to 
the initial-state electron but also to the final-state 
untagged electron. We use an exponentiation tech¬ 
nique [55] for the photon emission based on the pa¬ 
rameter 77 = (2a/n)(log(Qi/m 2 ) — 1) and the prob¬ 
ability density for the photon energy distribution, 
dP(rk) oc r/T dr/- , where 


Tk = 


jp* 

^ISR 


E, 


(14) 


beam 


As an approximation, the photon is always emitted 
along the incident electron direction on the tagged 
side. We limit the fractional energy of radiation to 
below r™ ax = 0.25 in the MC generation. 

In this configuration, the correction factor 1 + S 
to the tree-level cross section is close to unity |23| . 
Most of the events with large 77, typically 77 > 
0.1, are rejected by the selection criterion that uses 
-Eratio (the definition and criterion being described 
in Sec. |IV A I by requiring energy-momentum con¬ 
servation between the initial- and final-state particle 
systems without radiation. This effect is accounted 
for as a loss of efficiency for events with 77 < r™ ax . 

We generate events with a virtuality of the tagged- 
side photon Q 2 distributed with a constant form fac¬ 
tor over its continuous range Q 2 > 3.0 GeV 2 . The 
Q 2 value of each event is modified by ISR. We use the 


In this section, we describe the event selection 
and present some raw distributions to compare with 
those from MC. 


A. Selection criteria for signal candidate events 

A signal event contains an energetic electron and 
four photons. The kinematical variables are cal¬ 
culated in the laboratory system unless otherwise 
noted; those in the e + e - or 7*7 c.m. frame are 
identified with an asterisk in this section. We re¬ 
quire exactly one track that satisfies p t > 0.5 GeV/c, 
—0.8660 < cosf? < 0.9563, dr < 1 cm, and \dz\ < 
5 cm. There must be no other tracks that satisfy 
p t > 0.1 GeV/c, dr < 5 cm, and \dz\ < 5 cm 
in the above angular range. Here, p t is the trans¬ 
verse momentum in the laboratory frame with re¬ 
spect to the positron beam axis, 6 is the polar an¬ 
gle of the momentum direction with respect to the 
z axis, and (dr, dz) are the cylindrical coordinates 
of the point of closest approach of the track to the 
beam axis. We also require one or more neutral clus¬ 
ters in the ECL, whose energy sum is greater than 
0.5 GeV, as a pre-selection criterion for the exper¬ 
imental samples. These conditions are efficient in 
selecting a signal process within the kinematical re¬ 
gions of e + e _ —> e(e) 7 r° 7 r° in which one electron es¬ 
capes detection at small forward angles. 

For electron ID, we require E/p > 0.8 for the can¬ 
didate electron track. The absolute value of the 
momentum of the electron must be greater than 
1.0 GeV/c, where the electron energy is corrected for 











photon radiation or bremsstrahlung in the following 
way. In a 3° cone around the track, we collect all 
photons in the range 0.1 GeV< E 1 < p e c/3, where 
p e is the measured absolute momentum of the elec¬ 
tron track. The absolute momentum of the electron 
is replaced by p e + E E 7 . The cosine of the polar 
angle for the electron ( 0 e ) must be within the range 
—0.6235 < cos 9 e < 0.9481, which is the sensitive 
region for the HiE and Clst4 triggers. 

We search for a 7 r° candidate reconstructed from 
a photon pair with each photon having an energy 
above 0.1 GeV and a polar angle in the range 
—0.6235 < cos $ 7 < 0.9481. We constrain the po¬ 
lar angle of the photons from at least one 7r° in 
the sensitive region of the ECL triggers by the lat¬ 
ter condition, in order to reduce the systematic un¬ 
certainty of the trigger efficiency. The two-photon 
invariant mass is required to satisfy the criterion 
0.115 GeV/c 2 < M 77 < 0.150 GeV/c 2 . This tt° can¬ 
didate is referred to as “7rl”. 

In parallel, we search for 7 r° candidates with the 
7r°-mass-constrained fit among all pairwise combina¬ 
tions of photons in the entire ECL region. We select 
only combinations of photon pairs whose goodness 
of fit satisfies y 2 < 16. If ttI is also selected by the 
mass-constrained fit, we replace the four-momentum 
of 7rl by that of the result of the fit. 

We require that only one more pion be found 
among the n° candidates from the mass-constrained 
fit which does not share any photons with 7 rl. We 
refer to the second 7r° as u 7t2”. 

If there are two or more possible assignments of 7 rl 
and 7t2, we choose the one with the highest-energy 
photon to construct wl. If there are still two or more 
combinations that share the highest-energy photon 
in 7r 1, we choose the one in which the other photon 
in 7 r 1 has the higher energy. 

Figure [3] shows the 77 invariant-mass distribu¬ 
tion when constructing 7rl with a looser criterion of 
the two-photon invariant mass. The experimental 
data are compared with the distribution from the 
signal-MC sample using fits described below. The 
signal-MC distribution is fitted by the sum of a Crys¬ 
tal Ball function [25J and a linear function (signal 
component). Then, the experimental distribution 
is fitted by the sum of the signal component with 
the determined shape parameters and an additional 
linear function (background component), where the 
normalization and horizontal position for the signal 
component are allowed to float (resulting in a shift 
of the peak position by —0.8 MeV/c 2 ). The sole pur¬ 
pose of this fit is to compare the figures, which in¬ 
dicates a reasonable agreement and provides an es¬ 
timate of the background contamination. 
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FIG. 3: (a) The experimental distributions for the 77 

invariant mass forming 7 rl. The arrows indicate the se¬ 
lection range, (b) The distribution from the signal MC 
for the same variable. Statistics of the MC figure are 
arbitrary. The fit for the comparison is shown by a solid 
curve in each of the distributions (see the text). The 
background component in the experimental data from 
the fit is shown by the dashed line. 




We apply additional selection criteria for 7 rl and 
7 t 2 to reduce contamination from low-energy back¬ 
ground photons. The energy asymmetry for the two 
daughters ( 7 I and 72 ) of either pion, defined as 


_ | Eji — E 1 2 1 
E 7l + E 7 2 


(15) 


must satisfy E asym < 0.8. We require that the 
7 T° energies and transverse momenta satisfy E^i > 
0.4 GeV, E n2 > 0.3 GeV, E nl + E„ 2 > 1 GeV, 
Pt,ni > 0.15 GeV/c, and pt^ 2 > 0.15 GeV/c, respec¬ 
tively. We require that the polar angle of 7 t 2 in the 
laboratory frame satisfy —0.8660 < cos0 < 0.9563. 

We reject events with a back-to-back configura¬ 
tion of an electron and 7rl in the e + e~ c.m. frame, 
to suppress Bhabha events in which a track is not 
reconstructed; we require £*(e, 7 rl) < 177°, where 
C*(e, 7 rl) is the opening angle between the electron 
and the 7rl system. 

We require the tagged lepton to have the correct 
charge sign (“right-sign”) with respect to the beam 
from which it originates in the e + e - c.m. frame: 


9tagx(p: ie +p! >wl +p: iW2 )<o, (i6) 

where gtag is the tagged lepton charge. 

We apply a kinematical selection of 0.85 < 
-E-ratio <1.1, where £j- a ti 0 is defined as 


E r a fin — 


E 


♦ measured 


♦expected 


E 


(17) 


and £™ e 0 asured (£ , *® X o ected ) is the e + e c.m. en- 
ergy of the 7r°7r° system measured directly (expected 
by kinematics without radiation). This requirement 
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is motivated by a three-body kinematical calcula¬ 
tion for e + e _ — > e(e)R that is to be followed by 
R — > 7r°7T 0 , where R need not be a physical res¬ 
onance because this is a kinematical calculation. 
We impose a four-momentum conservation condi¬ 
tion Pinitiai(e + e“) = Pfi n ai(e(e)i?) wherein the direc¬ 
tion of the R momentum is taken to be parallel to 
that of the observed 7r°7r° system in the e + e _ c.m. 
frame. The expected energy of the R(= sys¬ 

tem, £X P 0 ected , is obtained by assigning the mea¬ 
sured 7r°7r° invariant mass to the R system. 

In Fig. [4j we show a comparison of E ratio be¬ 
tween the data and signal MC, where we observe a 
sharp peak corresponding to the signal process in the 
data that is consistent with the MC. The distribu¬ 
tion for “wrong-sign” events that have the opposite 
lepton charge to Eq. (16) is also shown; here, only 
a small peak is seen for the wrong-sign events near 
E ra tio = 1- This means that the backgrounds from 
e + e~ annihilation events, where the charge asym¬ 
metry of the tracks is not expected, are negligibly 
small. Meanwhile, there are significant right-sign 
backgrounds with a small E rat i 0 - We discuss such 
events in Sec. IV El 

We require transverse momentum balance in the 
e + e _ c.m. frame, |£pj| < 0.2 GeV/c, where 


m\ = \Pt,e+Pt,nl+Pi,2\- ( 18 ) 

We show the distribution for (referred to as 

ll p t balance”) in a wider range than the signal re¬ 
gion in three different W ranges in Fig. [5] where the 
samples after the three-body kinematics condition 
by the E lat io selection criterion applied are shown. 
The experimental distributions are compared with 
those of the signal-MC events. The signal peak near 
|Epj | = 0 is a little wider in the data than in the MC. 
This is partially due to the backgrounds in the data 
and partially due to non-inclusion of the finite-angle 
initial state radiation (ISR) on the tag side or the 
non-tag side. However, it is expected that the trans¬ 
verse momentum of the observed system, lEpJ'l, is 
dominated by those of the colliding virtual photons, 
which are taken into account in the MC simulation 
because the Q 2 of the electron after the ISR emis¬ 
sion is expected to be smaller than Q 2 of the colliding 
photon emitted by the electron. This is supported 
by our previous study of the |Epj| distribution for 
the 7*7 —>■ 7 r° process [ 10 ] . 

We find that events with .E rat ; 0 < 0.7 do not peak 
at | = 0. These features show that some charged 
or neutral hadrons or photons from 7r° decay es¬ 
cape detection in these events. They are consid¬ 
ered to come from multi-hadron production in two- 
photon processes or virtual pseudo-Compton scatter- 



E ratio 


x 10 


CM 

o 

o 

(co 

"E 

0 

> 

0 


0 

.Q 

E 

3 



FIG. 4: The distributions for Emtio for (a) the experi¬ 
mental data and (b) the signal-MC sample. The cross 
plots in (a) show the distribution for the wrong-sign 
events, and the arrows indicate the signal region. Statis¬ 
tics of the MC figure are arbitrary. 


ing (e + e _ —* e(e) 7 *, 7 * —► hadrons) according to the 
observed asymmetry in the correct- and incorrect- 
charge sign events, with a leakage of these back¬ 
ground components into the E rat i 0 signal region. 

The E ra tio and | Ep^ | distributions are discussed in 
more detail for background estimation, in Sec. [V] 


B. Assignment of kinematical variables for a 
signal event 

We assign four kinematical variables Q 2 , W, 
| cos 0*|, and tp* to each signal candidate event. The 
angles | cos 9* | and p* are defined in the 7*7 c.m. 
frame, where the direction of 7 * has cos 9* = 1 and 
the azimuthal direction of the recoiling electron de¬ 
fines <p* = 0 (Fig. [2|. Note that only these two vari¬ 
ables are defined in the 7*7 c.m. frame. In contrast, 
all the other variables with an asterisk appearing be- 
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FIG. 5: (a,b,c) The distributions for the pt balance for 
the data samples in three different W ranges indicated 
in each panel. The arrows show the selection region. 
(d,e,f) The corresponding distributions from the signal 
MC. Statistics of the MC figures are arbitrary. 



FIG. 6: The experimental distributions for W of the sig¬ 
nal candidates for 3 GeV 2 (5 GeV 2 ) < Q 2 < 30 GeV 2 
for the e-tag (p-tag) samples. Backgrounds are not sub¬ 
tracted. The cross and asterisk plots are for the e-tag 
and p-tag samples, respectively. 


C. Comparison of the experimental candidates 
with the signal-MC events 


low are defined in the e + e c.m. frame. 

The negative of the invariant mass squared, Q 2 , 
of the virtual incident photon is calculated using the 
measured four-momentum of the detected electron 
( p e ) from 

Qrec = (Pbeani Pe ) 

= 2E^ e&m E *(1 + q tag cos 0*), (19) 

where pbeam is the nominal four-momentum of the 
beam particle with the same charge as the detected 
electron; the right-hand side is given by the beam 
energy E £ eam and the observables of the tagged elec¬ 
tron in the e + e _ c.m. frame. We do not apply a cor¬ 
rection for initial-state radiation (ISR) on an event- 
by-event basis; instead, this effect is taken into ac¬ 
count in the correction for the differential cross sec¬ 
tion, as mentioned in Sec. |VI| 

The c.m. energy of the incident 7*7 collision, W, 
is the invariant mass of the final-state 7r°7r° system. 
The pion scattering angle 9* is defined in the 7*7 
c.m. frame as an angle between the virtual photon 
and that of one of the produced pions. In case an 
ISR photon is emitted in the tagged-electron side, 
the direction of the virtual photon is slightly misre- 
constructed and induces an error in cos 0*. However, 
the change due to this effect is typically ±0.01 in 
cos 9*. This is smaller than the angular resolution, 
typically cr cos g* = 0 . 02 , and the bin width, 0 . 2 , and 
thus we neglect the effect. 

There is a two-fold symmetry for the two sides of 
the plane, since (<£>*) = ~ ( P*)i so this 

angle is limited to the range 0 < ip* < 7 r. 


In this subsection, we show various distributions 
of the selected signal candidates. Backgrounds are 
not subtracted in the experimental data. Data dis¬ 
tributions are compared with the signal MC, where 
a uniform angular distribution and a representative 
Q 2 dependence are assumed. As most Q 2 and angu¬ 
lar dependence in these figures arise from kinemat¬ 
ics, such comparisons are meaningful, but no per¬ 
fect agreement between data and MC should be ex¬ 
pected. 

The experimental W distribution is shown in 
Fig. [ 6 ] for W < 2.5 GeV. For comparison, the corre¬ 
sponding distributions from the signal MC are also 
shown in the following figures. In Figs. [7] [9] all 

events within Q\ < Q 2 < 30 GeV 2 and 0.5 GeV 
< W < 2.1 GeV are integrated, where Q 2 L = 3 GeV 2 
(= 5 GeV 2 ) for the e-tag (p-tag) sample. We do 
not use the data below W < 0.5 GeV because the 
signal efficiency and the signal-to-background ratio 
decrease steeply below that energy. No large dis¬ 
crepancy between the data and signal MC is seen in 
all the figures. This implies that the MC is a faith¬ 
ful representation of reality and the backgrounds in 
the experimental data are not very large. The dif¬ 
ference in the pion energy distributions reflects the 
difference in its angular distribution in the 7*7 c.m. 
frame between the data and signal MC, with the uni¬ 
form angular distribution for the latter. 


V. BACKGROUND ESTIMATION 

We consider several sources of possible background 
processes that could be misidentified as the signal 
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FIG. 7: (a,b,c) The Q 2 distributions for the data samples 
in three different W ranges indicated above each panel. 
The cross and asterisk plots are for the e-tag and p-tag 
samples, respectively. (d,e,f) The corresponding distri¬ 
butions from the signal MC, where the solid and dashed 
histograms are for the e-tag and p-tag samples, respec¬ 
tively. Statistics of the MC figures are arbitrary, but the 
scale is common for the e- and p-tags so their ratio could 
be compared between MC and data. 
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A. Single pion production process 

Backgrounds from the single pion production, 

e+ e - e ( e ) 7r 0 ) 

with one fake 7r° that is wrongly re¬ 
constructed, are estimated by MC simulation of the 
process. From our previous measurement m , we 
know the cross section of this process with a suffi¬ 
cient accuracy for this purpose. From the MC study, 
the contamination is estimated to be less than one 
event in the whole sample of candidate events (about 
3700 events) and is enhanced in the forward angular 
region in the c.m. frame, as shown in Fig.JlOjc). As 
the estimated number of events is small, we neglect 
this background source. The Q 2 dependence of this 
background source is expected to be similar to that 
of the signal process. 


B. Radiative Bhabha process 

The radiative Bhabha process with the virtual 
Compton scattering topology, e + e - —> e(e)j, which 
has a relatively large cross section, has been proved 
to give only a small contribution to the present mea¬ 
surement. This is verified in the 7rl and n2 mass 
distributions (Fig. [8]), where the pion peaks in the 


FIG. 8 : (a,b,c,g,h,i) The distributions for the 77 invari¬ 
ant mass to construct 7rl (a,b,c) and 7t2 (g,h,i) for the 
data samples in three different W ranges indicated above 
each panel. (d,e,f,j,k,l) The corresponding distributions 
from the signal MC for 71 T (d,e,f) and 7r2 (j,k,l). The 
legend and comments are the same as those in Fig. [7] 

experimental data with a high purity compared to 
the signal-MC samples. This background must be 
less than 5% according to the shape of the m n 1,2 
distribution. However, we must note that this back¬ 
ground forms a broad enhancement near the pion 
mass when the high-energy photon is converted to 
an e + e _ pair in front of the ECL. The measured Q 2 
values for the virtual Compton scattering tend to 
populate the region of higher values than in the sig¬ 
nal process, and the background should be relatively 
large in the liigh-Q 2 region. 

The estimation of the contamination from this 
background process using a background-MC sample 
is very difficult because of the large MC statistics 
needed due to the large cross section of the process 
and the significant suppression in the event selection. 
Thus, we estimate the contamination using distribu¬ 
tions of the experimental signal candidates. 

In presence of both of the described background 
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FIG. 9: (a,b,c) The distributions for the laboratory en¬ 
ergy of 7r 1 and 7 t 2 for the data samples in three different 
W ranges indicated above each panel. (d,e,f) The cor¬ 
responding distributions from the signal MC. (g,h,i) and 
(j,k,l) are the distributions for the laboratory angle of 7rl 
and 7 t 2 for the data samples in the indicated W ranges 
for the experimental samples and MC samples, respec¬ 
tively. The legend and comments are the same as those 
in Fig. 0 


processes, single pion production and the radiative 
Bhabha process, at least one of the two 7r°s recon¬ 
structed is not a true n°. A fake ir° is wrongly 
reconstructed from the beam background or elec¬ 
tronic noise. If such a background n° had large re¬ 
constructed energy, the event would have been re¬ 
jected by the p t -balance selection criterion. There¬ 
fore, these fake pions remain in the signal event sam¬ 
ple only if the noise has small reconstructed energy. 
In addition, such noise photons populate the lowest 
energy region just above the energy thresholds for a 
photon or a neutral pion for the selection. Conse¬ 
quently, the c.m. scattering angle in the two-photon 
system tends to be reconstructed in the forward re¬ 
gions because of the unbalanced laboratory energies 
between the two pions and the Lorentz boost of the 


two-pion system along the beam axis, just as ex¬ 
pected in the single-7r° production case. However, 
as seen in Fig. [lO|a), no visible enhancement is seen 
even in the most forward region | cos0*| >0.9, where 
the number of the signal events is small according to 
the small signal acceptance (Fig. [Hb)) in the an¬ 
gular distribution for the entire experimental signal 
candidates. Figure |Tl|a) shows the 7rl and 7 t 2 mass 
distribution for the events with | cos d* \ > 0.9 in com¬ 
parison with the signal-MC expectation and the ex¬ 
perimental data for | cos 9* \ < 0.9 (Fig. 0b) and 
(c), respectively). It is difficult to estimate the non- 
7 r° backgrounds quantitatively from these figures. 

In our previous study m, we found that the 
backgrounds from virtual Compton scattering have 
a high-energy neutral pion from misreconstruction 
of an e + e~ conversion of the single photon, which 
simulates two photons with a high probability. Its 
characteristic feature is a small polar angle differ¬ 
ence AO due to the effect of the magnetic field. It 
is found experimentally that the backgrounds con¬ 
centrate strongly at Ad E^ < 0.05 rad-GeV [TU]. 
Figure 12 shows the distribution of A 0 E^ for nl. 
We estimate the yield of the backgrounds from the 
virtual Compton process to be around 10 events, 
which are seen as an enhancement in the region 
Ad E^ < 0.05 rad-GeV, out of the 67-event subsam¬ 
ple with | cos d* | > 0.9 and 0.5 GeV< W < 2.1 GeV. 
We thus estimate that the background fraction is 
about 15% for | cos 0*| > 0.9, although its uncer¬ 
tainty is large. We do not find any sign of back¬ 
ground from this source in the region | cosd*| < 0.9. 


C. 7T°7 production process 

Backgrounds from the 7r°7 production, e + e - —> 
e(e)7r°7, are dominated by the virtual pseudo- 
Compton scattering process of w-meson production, 
e + e _ —> e(e)ui. We know the cross section of this 
process with sufficient accuracy from our previous 
measurement m to estimate this contamination 
from the background-MC study, as described below. 

The W, Q 2 , and |cos0*| dependences are based 
on the MC. Figure [13] shows the distributions of 
the background from the process e + e~ —> e(e)u>, 
to —>■ 7r°7 in the Q 2 -W, and |cos(9*|-IV directions. 
This background is enhanced in the W region near 
the uj mass (W = 0.7 - 0.9 GeV) in the high¬ 
ly 2 region ( Q 2 > 12 GeV 2 ). The total number of 
events contaminating the signal process is estimated 
to be about 9 events, of which only 6 events for 
Q 2 > 12 GeV 2 are expected to have a non-negligible 
effect compared to the signal yield. The | cos 9* | dis- 
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tribution is approximately flat so the behavior of this 
background is very different from the former back¬ 
ground sources where both photons constituting one 
pion originated from noise. We find no visible en¬ 
hancement in the 0.7 - 0.9 GeV region in the ob¬ 
served W distribution, as expected from the small 
contamination of the background process. 


rf —> 7r°7r°7r 0 is separately estimated using the rf- 
TFF El. the branching fraction, and generated 
background-MC samples for the process. We find 
that the expected contribution for W > 0.5 GeV is 
about 0.3 events. Thus, we conclude that the contri¬ 
bution of this background is negligibly small. 


D. Three- x° production process 


The three-pion production process, 7*7 
7r 0 7r°7r°, is wrongly selected as the signal candidate if 
one non-energetic pion escapes detection. As three- 
pion production in single-tag two-photon collisions 
has not been measured to date, we estimate the con¬ 
tamination by referring to the corresponding zero- 
tag measurement. 

The 7r + 7r _ 7r° production in the zero-tag process 
has been measured by the L3 collaboration [ 22 ] 
and shows that the process is dominated by the 
02(1320) and a structure near 1.7 GeV (02(1700) and 
potential 7^(1670) production decaying to p°ir° or 
/2(1270)7 t°). However, the three-neutral-pion pro¬ 
duction is strongly suppressed because p° does not 
decay to 7r°7r° and 02(1700) —* /2(1270)7 t° has not 
been observed definitively. The L3 measurement 
only provides an upper limit for 7^(1670) production 
from the two-photon process. A measurement of a 
finite value for 77 —► 7t 2 (1670) —► 7r 0 7r 0 7r° production 
is reported by the Crystal Ball experiment m but 
it is not consistent with the L3 upper limit. 

We estimate that 7r 0 7r°7r° production is about 4% 
of 7r°7r° production in the zero-tag two-photon pro¬ 
cess 1211122] in the e + e _ -based cross section, accord¬ 
ing to the measurement of L3 for 77 —> 02 (1700) —► 
/2(1270)7t° [26| . We confirm this estimate by our 
count of zero-tag 77 —► 7r°7r n 7r 0 candidate events, 
requiring three neutral pions in an event, from the 
Belle data samples. Our estimate of the cross section 
is somewhat smaller than that from the above Crys¬ 
tal Ball measurement E2, but is consistent with it 
within about a factor of two. 

We assume that the cross section ratio of the 
two processes, 7r 0 7r°7r° to 7r°7r°, is the same between 
the zero-tag and single-tag processes. Taking into 
account the selection efficiency for the background 
events and the background subtraction using the p t - 
balance distribution (applied in Sec. VIB), the esti¬ 
mated contamination is less than 0.5% of the signal 
yield, where the 7r 0 7r°7r° background does not peak 
near |Ep(| = 0, like for the other non-exclusive pro¬ 
cesses described in the following section. 

The contribution from the process 7*7 —» 


E. Other non-exclusive processes 

The other non-exclusive background processes, 
e + e _ —► e(e)7r°7r°A', where X denotes mul¬ 
tiple hadrons, are in general subdivided into 
two-photon (C-even) and virtual pseudo-Compton 
(bremsstrahlung, C'-odd) processes, but they inter¬ 
fere with each other if the same X is allowed for 
both processes. The majority of such background 
events populate the small E rat i 0 region, e.g., less than 
0.7. This feature is distinct from the aforementioned 
background processes that can populate the region 
near E rat i 0 = 1. These backgrounds also do not peak 
near | = 0. 

From Fig. [4] we expect that the low-E rat i 0 compo¬ 
nent could leak into the signal region. We estimate 
the relative ratio of the non-exclusive backgrounds 
to the signal yield by counting the number of events 
in the subregion of the signal region (0.85 < E rat i 0 < 
0.925 n 0.1 GeV/c < |Ep*| < 0.2 GeV/c) where 
the background component would be relatively large. 
We assume that the backgrounds are distributed as 
a linear function in 0 < |Hpt| < 0.2 GeV/c and 
in 0.85 < E ratio < 1, vanishing at |£p*| = 0 and 
E ra tio = 1, respectively, factorized for the two direc¬ 
tions. The fraction of the backgrounds falling in the 
above subregion is calculated to be 9/16 of that in 
the entire signal region, according to the distribu¬ 
tion. We also estimate the fraction of signal events 
coming into the same region using the signal-MC 
sample. We thus determine the expected number 
of the background events from this information. 

In addition, we use the W and Q 2 dependence of 
the E rat io sideband events, which are extracted from 
the experimental data in 0.7 < E rat i 0 < 0.8 as an¬ 
cillary information for the kinematical regions with 
low statistics. The normalization of the background 
is determined in the first method in a high-statistics 
region; we extrapolate it to different W and Q 2 re¬ 
gions with lower statistics assuming the dependence 
observed in the E rat i 0 sideband region. 

We show the E rat i 0 and | | distributions in the 

signal region of the experimental and signal-MC sam¬ 
ples in Figs. [14] and [15] respectively. The dashed 
lines in the panels for the experimental data show 
the estimated non-exclusive background. 
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FIG. 10: The |cos0*| distributions for (a) the exper¬ 
imental signal candidates, (b) the signal MC with the 
isotropic generation for | cos0*|, and (c) the background- 
MC events for the single-pion production process, where 
the yield is scaled to the expected contamination in the 
set of all signal candidates. The error bars are statistical 
from the MC samples but are not proportional to square 
root of the number of events because each event in the 
MC sample has a variable weight to reproduce the ex¬ 
perimental Q 2 dependence. For (b), the normalization is 
arbitrary. 



mn 2 (GeV/c 2 ) 


FIG. 11: Scatter plot for the 77 invariant masses 
for 7r 1 and n2 for (a) the experimental signal candi¬ 
dates for | cos 0*| > 0.9, (b) the signal-MC samples for 
| cos 0*| > 0.9, and (c) the experimental signal candidates 
for | cos 0*| < 0.9. 


Exp. |cos 0*|>O.9 


Exp. |cos 0*|<O.9 


> 

CD 

CD 



Signal MC |cos 0*|>O.9 


Signal MC |cos 0*|<0.9 



FIG. 12: The distributions of A0E 77 for (a) the exper¬ 
imental signal candidates in |cos0*| > 0.9 and (b) in 
| cos 0*| < 0.9. Similar plots for (c) the signal-MC sam¬ 
ples in | cos0*| > 0.9 and (d) | cos0*| < 0.9. The statis¬ 
tics of the MC figures are arbitrary. 
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FIG. 13: The distribution of background events esti¬ 
mated from the MC for the process e + e~ —> e(e)ui, 
uo —> 7r°7 in the Q 2 versus W and | cos 0* | versus W scat¬ 
ter plots. About three events in each plot correspond 
to a contamination of one event in the present signal- 
candidate sample. 


VI. DERIVATION OF THE DIFFERENTIAL 
CROSS SECTION 


The details of the background subtraction are de¬ 
scribed in Sec. IVIBl 


We first define and evaluate the e + e _ -incident- 
based cross section separately for the p-tag and e- 
tag samples. Then we derive the differential cross 
section of the process 7*7 —> 

The e + e _ -incident-based differential cross section 
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Exp. Signal MC Exp. Signal MC 



FIG. 14: Comparison between the experimental (left) 
and signal-MC (right) distributions for i? ra tio of events 
in the signal regions. Each row corresponds to the 
same IT and Q 2 regions indicated in the right panel. 
The dashed lines show the estimated non-exclusive back¬ 
ground, which is assumed to distribute linearly with the 
horizontal variable. Statistics of the MC figures are ar¬ 
bitrary. 


FIG. 15: Comparison between the experimental (left) 
and signal-MC (right) distributions for the pt balance of 
events in the signal regions. Each row corresponds to 
the same IT and Q 2 regions indicated in the right panel. 
The dashed lines show the estimated non-exclusive back¬ 
ground, which is assumed to distribute linearly with the 
horizontal variable. Statistics of the MC figures are ar¬ 
bitrary. 


is written as 


( d' <j ee \ 

ydWd I cos 9*\dQ 2 Is—tag 

_tx-taglw.icoseTQ 2 )_ 

e(, (W,| cos 0* |,Q 2 )AWA| cos 0* | A Q 2 f CdtB 2 


, ( 20 ) 


where the yield Y and the uncorrected efficiency 
obtained by the signal MC e' are separately eval¬ 
uated for p-tag and e-tag, for a consistency check. 
The measurement ranges of W, | cos 9* |, and Q 2 
and the corresponding bin widths AIT , A|cos0*|, 
and A Q 2 are summarized in Table III where the 
differential cross section was first calculated with 
AIT = 0.05 GeV and A| cos0*| = 0.1, and then two 
adjacent IT and |cos0*| bins are combined (for IT 
only above IT > 1.1 GeV) with an arithmetic mean. 
Here, J Cdt is the integrated luminosity of 759 fb -1 
and B 2 = 0.9766 is the square of the decay branching 
fraction B(tt° — > 77). 

We take into account the difference of the beam 
energies in evaluating s'. Since the efficiency and 
the luminosity function depend on the beam energy, 
we construct the corrected efficiency using the ratio 


of the products of the efficiency, luminosity function, 
and the integrated luminosity for two cases: that the 
products are combined for the different beam ener¬ 
gies in the experiment and that all the experiment 
would be done at the sole T(4S') energy with the 
same total integrated luminosity as in the experi¬ 
ment. Thus, the e + e _ -based cross section measured 
for the energy of T(45'), 10.58 GeV, is obtained. 

After confirming the consistency between the p- 
and e-tag measurements to ensure validity of the 
efficiency corrections (described in Sec. VIA I, we 
combine their yields and efficiencies using the for¬ 
mula which builds in the equality of the efficiency- 
corrected yields for both measurements, 


_ d 3 a ee _ _ 

dWd\ cos 6* | dQ 2 

Y(W, | cos0*|,Q 2 )(1-6(W,| cos 6* \ ,Q 2 )) 
e'{W,\ cos0*|,Q 2 )AWA| cos0*|AQ 2 / Cdt ’ 


( 21 ) 


where Y = T p _ tag + T e -t ag , e[ = (£p-t ag + £ e-t ag )/ 2 . 
and b is the background fraction combined for p- and 
e-tags, which is subtracted here. For the region Q 2 = 
3-5 GeV 2 , we do not use the p-tag data because its 
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TABLE III: The measurement range and bin widths for 
three-dimensional variables (IF, | cosf?*|, Q 2 ). 


Variable 

Measurement 

range 

Bin width 

Unit 

Number 

of bins 

IF 

0.5 - 1.1 

0.05 

GeV 

12 


1.1 - 2.1 

0.1 


10 

| cos 9* | 

0.0 ~ 1.0 

0.2 


5 

Q 2 

3.0 - 5.0 (e-tag only) 

1.0 

GeV 2 

2 


5.0 - 6.0 

1.0 


1 


6.0 - 12.0 

2.0 


3 


12.0 15.0 

3.0 


1 


15.0 - 20.0 

5.0 


1 


20.0 - 30.0 

10.0 


1 


statistical accuracy is much worse than for the e-tag 
sample. There, as a result, the cross-section value in 
the e-tag measurement is simply doubled. 

Finally, the e + e _ -incident-based differential cross 
section is converted to that based on 7* 7-incident 
by dividing by the single-tag two-photon luminosity 
function d 2 L 1 * 1 /dWdQ 2 , which is a function of W 
and Q 2 . We use the relation 


da. 


7 7 


d | cos 6 *| 


dWd \cosS W ^_ iil+5){e/s , )£ , 


( 22 ) 


The factors 6, e, and / correspond to the radiative 
correction, efficiency corrected for the p* dependence 
of the differential cross section, and the unfolding 
effect that accounts for migrations between the dif¬ 
ferent Q 2 bins, respectively, which are explained in 
more detail in the subsections below. 


A. Efficiency plots and consistency check of 
the p-tag and e-tag measurements 

Figure [lB] shows the trigger efficiencies obtained 
from the signal-MC samples and trigger simulator. 
The trigger efficiency is defined for events within the 
selection criteria. The W dependence of the trigger 
efficiency is mild in the measurement region. The 
dip-bump structure seen in the Q 2 dependence for 
the p-tag efficiency is an artifact of the Bhabha-veto 
logic in the HiE trigger. 

Figure [T7] shows the efficiencies in which all the se¬ 
lection and trigger conditions are taken into account. 
They are provided as a function of W and | cos 9* | for 
the selected Q 2 bins of the p- or e-tag samples. These 


(a) W = 0.6 GeV 



f (b) W = 1.0 GeV (c) W = 1.3 GeV 



(d) W = 0.6 GeV 



(e) W= 1.0 GeV 



(f) W = 1.3 GeV 



FIG. 16: Trigger efficiency estimated by the signal-MC 
samples and trigger simulator. The solid (black) and 
dashed (red) histograms are for e-tag and p-tag events, 
respectively. The results are for the three W points in¬ 
dicated above each panel: (a,b,c) the Q 2 dependence for 
the isotropically generated 7r° pairs in the 7*7 c.m. sys¬ 
tem; (d,e,f) the scattering- (polar-) angle dependence 
of 7 T° in the 7*7 c.m. system. The Q 2 range 3 GeV ' 2 
(5 GeV 2 ) < Q 2 < 30 GeV 2 is integrated for the e-tag 
(p-tag) plot. 


efficiencies are obtained from the signal-MC events, 
which are generated assuming an isotropic angular 
distribution of the pions in the 7*7 c.m. frame. Ef¬ 
ficiency corrections for the <p * dependence are not 
taken into account in these figures. 

Our accelerator and detector systems are asym¬ 
metric between the positron and electron incident 
directions and energies, and separate measurements 
of the p-tag and e-tag samples provide a good valida¬ 
tion check for various systematic effects of the trig¬ 
ger, detector acceptance, and selection conditions. 
Figures [l8|and[T9]compare the e + e _ -based cross sec¬ 
tion measured separately for the p- and e-tags. They 
are expected to show the same cross section accord¬ 
ing to the C symmetry if there is no systematic bias. 
The results from the two tag conditions are consis¬ 
tent within statistical errors. 


B. Background subtraction 

We have estimated the background yields in the 
signal samples using the background-process-MC 
samples for the single pion production, e + e _ —> 
e(e)j r°, and the 7r°7 production dominated by 
e + e _ —► e(e)cu, uj —1- 7r°7. These are described in 
Ref. f 10] , where their yields are normalized based on 
the observations. 
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FIG. 17: Efficiency estimated from the signal-MC samples passed through all the selection and trigger conditions. 
Each panel provides the W and | cos 9* | dependences of the efficiency for the four selected Q 2 bins for each of the p- 
or e-tags. The contour levels shown by different colors are not common across the panels. 


As mentioned in Sec.[V] the estimated total yields 
for the sum of p- and e-tags are less than one and nine 
events for the single pion and the 7r°7 production 
backgrounds, respectively. 

The single-pion background is small, and we do 
not subtract its contribution but rather include its 
effect in the systematic uncertainty due to subtrac¬ 
tion of the virtual Compton scattering. 

For the 7r°7 production, the background is esti¬ 
mated to be 2%, 6%, and 13% of the signal candi¬ 
date yields for the Q 2 bins of 12 - 15 GeV 2 , 15 - 
20 GeV 2 , and 20 - 30 GeV 2 , respectively, indepen¬ 
dently of W and | cos0*|, and is subtracted. For the 
bins at Q 2 < 12 GeV 2 , we neglect this background 
source. 

The details of the estimation of the background 
contributions from the virtual Compton process are 
discussed in Sec. |V B| We estimate that the signal- 
candidate yield in each |cos0*| > 0.9 bin contains 
(15 ± 15)% of the background from this process, 
and we derive the differential cross section for the 
| cos 0*| > 0.8 bin. We estimate systematic uncer¬ 
tainties conservatively because possible Q 2 depen¬ 
dence and the background from the single-pion pro¬ 
duction process with similar properties are neglected 
here, and the contamination could be sensitive to the 
noise conditions. The error size is estimated from 
A 6 E 11 and | cos 6* | distributions for several differ¬ 
ent conditions. The uncertainties of the background 


estimation for different kinematical regions are dis¬ 
cussed in Sec. IVI FI 

We neglect the contribution of the 7r°7r 0 7r° produc¬ 
tion process because it is estimated to be less than 
1% of the signal process (Sec. VDl. 

The contamination of the other non-exclusive 
background processes is estimated using the E la f IO 
and | variables as described in Sec. VE We es¬ 
timate the fraction of the background in each of the 
(W, Q 2 ) regions to be between 3% and 12%. This 
background fraction has no prominent | cos 9* | de¬ 
pendence and we neglect it. 

In subtraction of the background yields, we mul¬ 
tiply the observed yield by the expected ratio of the 
background to the observed yields in each bin, as 
represented by the (1 — 6) factor in Eq. (211. 


C. Efficiency corrections 

The efficiency after the integration over ip* is cal¬ 
culated using the signal MC, assuming a flat ip* dis¬ 
tribution. We correct the efficiency according to the 
actual non-uniformity in ip* observed in the data. 
This correction is necessary in case both the effi¬ 
ciency (e) and the differential cross section (da/dip*) 
have a ip* dependence. 

We partition the kinematical region of the mea¬ 
surement into three-dimensional (W, |cos0*|, Q 2 ) 
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FIG. 18: The W dependence of the e + e _ -based cross section in each Q 2 bin. The full |cos#*| range (0 to 1) is 
integrated. The closed circles (open circles) are for the e-tag (p-tag) measurements. The plotting location for each 
p-tag point is shifted slightly to improve the visibility of error bars. The efficiency corrections and background 
subtraction are applied. 


rectangular-prism cells defined by 0.5 GeV < W < 
1.1 GeV, 1.1 GeV < W < 1.6 GeV, and 1.6 GeV 
< W <2.1 GeV; five equal-width bins in cos 0*|; 
3 GeV 2 < Q 2 < 8 GeV 2 , 8 GeV 2 < Q 2 < 12 GeV 2 , 
and 12 GeV 2 < Q 2 < 30 GeV 2 . 

Due to limited statistics, we take such a coarse 
binning for the W and Q 2 directions, taking into ac¬ 
count qualitative changes of | cos 6* | dependence of 
the differential cross section. We estimate the effi¬ 
ciency correction factor, e/e' , in each of the cells. 

The signal-MC distribution after the selection be¬ 
haves as 1 Vmc(v*) oc and the experimen¬ 

tal distribution as -/Vexp(v?*) cx £ { l P*)P( l P*), where 
P(p*) oc da/dtp*. Then, the efficiency correction 
factor is calculated as 

£ = ^So e{v*)P{v*)dy* 

e' So Jo P(<fW { 


The ip* dependence of the efficiency and the p*- 
differentiated cross section are obtained by e(p*) oc 
Nmc(</?*) and P(p*) oc N EX p{<P*)/N MC {p*), respec¬ 
tively. We then expand each function as a Fourier 
series, 

-Vmc( < P*) = -4(1 + c cos p* + dcos2p* + ...) (24) 

and 

N exp Sp* ) /-Vmc (^*) = B{ 1 + acosp* + 6cos2v?*), 

(25) 

where the coefficients are determined by fitting. 
There are no sine terms because we expect symmetry 
between p* and — p* for these functions, or we can 
regard them as the sum of the functions in the posi¬ 
tive and negative p* regions, g(p*) = f( t P*)+f(—'P*) 
(0 < p* <7r). They result in 


2 


(26) 
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FIG. 19: The | cos#*| dependence of the e + e _ -based cross section in each Q 2 bin. The dependence for three different 
W regions is shown in above each panel, and the W regions are 0.5 - 1.1 GeV, 1.1 1.6 GeV, 1.6 - 2.1 GeV from 

bottom to top within each panel, as indicated in the upper left panel. The closed circles (open circles) are for the 
e-tag (p-tag) measurements. The plotting location for each p-tag point is shifted slightly to improve the visibility of 
the error bars. The efficiency corrections and background subtraction are applied. 


This formula is independent of the normalizations 
of the Nmc(<P*) and -Vexp(</ 2*) functions. We ap¬ 
proximate A r Mc( < / 2 *) with a Fourier expansion up to 
cos Aip *; the coefficients of cos 3y>* and cos 4 <p* terms 
do not affect the e/e' result because P(<p*) is up to 


only cos 2 <p*, but the effect of the terms is signif¬ 
icant in the fit to determine the coefficients c and 
d. The ip* dependence of the efficiency and the fit 
for two W regions in the Q 2 range between 8 GeV 2 


and 12 GeV 2 is shown in Fig. 20 The experimen- 
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FIG. 20: The p* distributions for e(p*) oc Nmc{p*) 
for the W region 0.5 - 1.1 GeV (upper row) and 1.1 - 
1.6 GeV (lower row), and 8 GeV 2 < Q 2 < 12 GeV 2 . They 
are separately plotted for the five angular regions, whose 
central | cos 9* | value is indicated in each panel. The solid 
curve is the fit to the Fourier expansion described in the 
text. 


W = 0.5 -1.1 GeV W=1.1-1.6 GeV W= 1.6 -2.1 GeV 



|cos 0*1 


tal results for the ^A-differentiated cross section are 
discussed in Sec. IVII Cl 

We find that the ip* dependence of the differ¬ 
ential cross section, and thus the correction-factor 
value, change drastically between above and be¬ 
low | cos 0*| = 0.6 and that the value for each of 
| cos 9* | < 0.6 and | cos 9* | > 0.6 is almost constant 
for the cells with the same (W, Q 2 ). Thus, we take 
a weighted average of the correction factor for the 
three (two) bins of |cos0*| < 0.6 (|cos0*| > 0.6), 
in order to reduce the uncertainty of the correction. 
The obtained efficiency correction factor e/e', as well 
as the fit results for a and b , are plotted in Fig. [21] 
The correction factor ranges from 0.67 to 1.31 and is 
within 0.93 - 1.06 for 10 of the 18 regions. 

The trend of the ip * dependence is explained by the 
interference term(s) with the D 0 ( J = 2 and A = 0) 
component, which changes its sign at l/-\/3 (~ 0.577) 
according to the V 2 ° function (see Eq. ©)• 


FIG. 21: The | cos 6* \ dependence of the obtained param¬ 
eters of a, b, and the efficiency correction factor e/e '; a 
and b are the coefficients of cos p* and cos2<£* in P((p*) 
for each of the nine (W, Q 2 ) regions (see the text). Thick 
lines indicate the efficiency correction factor, which is 
applied to obtain the ip*-integrated cross section. 

TABLE IV: The nominal Q 2 value ( Q 2 ) for each Q 2 bin. 


Q 2 bin (GeV 2 ) 

Bin center (GeV 2 ) 

Q 2 (GeV 2 ) 

3-4 

3.5 

3.45 

4-5 

4.5 

4.46 

5-6 

5.5 

5.47 

6-8 

7.0 

6.89 

8-10 

9.0 

8.92 

10 - 12 

11.0 

10.93 

12 - 15 

13.5 

13.37 

15 - 20 

17.5 

17.23 

20 - 30 

25.0 

24.25 


D. Radiative correction and Q 2 unfolding 


We apply a correction of 2% (S = 0.02) for the 
cross section as the radiative correction. This is the 
same as that evaluated in the single pion produc¬ 
tion m ■ This correction depends very little on W, 
| cos 0*|, and Q 2 , and is treated as a constant. 

We define the nominal Q 2 for each bin with a finite 
bin width, Q 2 , using the formula 


da ee 

dQ 2 


m 



d(J ee 


( Q 2 )dQ 2 , 


(27) 


where A Q 2 is the bin width. We assume an approx¬ 
imate dependence of da/dQ 2 oc Q ~ 7 for the calcula¬ 
tion m, independent of W and | cos 0* |, and have 
omitted the notations of W and | cos0*| in Eq. (27). 


The Q 2 values for the Q 2 bins are listed in Table IV 
We use the luminosity function at this Q 2 point to 
obtain the 7*7-based cross section for each Q 2 bin. 
We also use the central value of the Q 2 bins to rep¬ 
resent the individual bins in tables and figures for 
convenience of description. 
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The Q 2 value measured for each event differs from 
the true Q 2 for two reasons: the resolution effect 
of the Q 2 determination and the reduction of the 
incident electron energy due to ISR. 

The signal yield is measured in bins of recon¬ 
structed Q 2 ec . For the true value of Q 2 , we use the 
corrected value Q 2 or , and we assign the number in 
Q 2 OI for each of the measurement bins. According to 
the ISR effect found in the signal-MC events, Q 2 ec 
is about 1% larger than Q 2 OI on average. We es¬ 
timate the displacement of events across bins using 
the signal-MC events in order to unfold the Q 2 dis¬ 
tributions folded by the displacement. 


TABLE V: Value of the eo parameter depending on Q 2 
for the two W regions. 


Q 2 bin (GeV 2 ) 

W region 
0.5 - 1.6 GeV 1.6 

- 2.1 GeV 

3-4 

0.82 

0.77 

4-5 

0.88 

0.84 

5-6 

0.90 

0.83 

6-8 

0.89 

0.83 

8-10 

0.88 

0.85 

10 - 12 

0.88 

0.85 

12 - 15 

0.86 

0.83 

15 - 20 

0.82 

0.80 

20 - 30 

0.76 

0.73 


We correct the measurement by the factor /) = 
Nji/ Nij, following the method applied to 
the Q 2 edge regions in our previous analysis for 
7 t°-TFF |T0], where N t j is the transfer matrix ob¬ 
tained from MC simulation with the adjusted in¬ 
dependence close to the observed dependence; here i 
(j) is the bin number for Q 2 ec (Q 2 or ). 

However, the correction for this Q 2 migration ef¬ 
fect is already partially included in the efficiency de¬ 
termination. We have defined the efficiency as the 
selected number of events in a Q 2 ec bin, as in the ex¬ 
perimental data, but divided by the number of gen¬ 
erated events in the Q 2 or bin. This is due to the 
difficulty to define Q 2 ec for generated events falling 
outside the detector acceptance. Thus, a migration 
effect is also introduced in the efficiency, but more 
weakly than in the experimental data, at half the size 
of the latter, due to the different Q 2 dependencies. 
By the consideration of this point, we correct the 
cross section for this effect by the factor (fi — 1)/2 
and assign a systematic uncertainty of the same size 
as the correction value. The correction is +5% for 
the lowest Q 2 bin and within ±3% for the other Q 2 
bins. We assign the systematic uncertainty for this 
correction with about half the size of the correction, 
2 %. 

We do not use the matrix inversion unfolding 
method m because the low statistics in the re¬ 
lated multidimensional bins would tend to enhance 
the statistical errors. In addition, a systematic bias 
could appear in the edge bins of the measured Q 2 
range by this method. 

The measured 7*7-based cross section cr 7 * 7 = 
cttt + eoO'LT = u is discussed and shown in figures 
in Sec. IVIII 


E. Effect of <tlt (helicity-1) component in the 
signal sample 


We estimate eo, the factor multiplying ctlt in 


Eq. (121, in each bin. We use the mean value of 
eo calculated by Eq. ([8]) for each selected event from 
the signal-MC samples in different kinematical re¬ 
gions. The value of eo has a weak dependence on 
Q 2 , W, | cos 6*\, and ip*. Since the W dependence, 
in W = 0.5 - 1.5 GeV and the |cos0*| dependence 
are small (within ±4%), we neglect their effect. It 
has some p* dependence (up to ±7%). 

We have compared the Q 2 dependence of the ex¬ 
perimental events with that of the signal-MC sam¬ 
ples and have confirmed their consistency. Similar 
calculations for are also performed. 

We tabulate the Q 2 dependence of eo in Table [v] 
for the two W regions. The expected values are eo = 
0.88 ± 0.06 and ei = 1.28 ± 0.07 for W = 1.1 - 
1.5 GeV for the Q 2 -integrated samples. The ranges 
in si and e 2 are due to the p* dependence. These 
values are used for partial-wave analyses performed 
with Q 2 -dependent and Q 2 -integrated cross sections 
in Sec. VIIC and [VII D[ respectively. 


For a quantitative study of the A = 1 component, 
we use the azimuthal-angle differentiated cross sec¬ 
tion d 2 cr/d\ cos 9*\dp*. It is derived as follows. 

The p* dependence of d 2 a/d\ cos 6*\dp* follows 
^exp^VAmc^*) in each W and | cos 0*1 bin in¬ 
tegrated in the Q 2 = 5 - 30 GeV 2 region. For this 
purpose, we use five W bins, 0.7 - 1.1 GeV, 1.1 - 
1.2 GeV, 1.2 - 1.3 GeV, 1.3 - 1.4 GeV, and 1.4 - 
1.5 GeV. 

The p* dependence is normalized to the arithmetic 
mean (da/d\ cos 9* |) of the differential cross section 
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over the corresponding (W, Q 2 ) bins. We do not use 
the results for Q 2 < 5 GeV 2 because poor statistics 
in that region would diminish the accuracy of the 
arithmetic mean. 

The normalization results in 

f 71 d 2 a ^ * _ / da N 

Jo d\ cos 9*\d<p* ^ ~ { d\cosd*\ } - 1 j 

The results for d 2 a/d\ cos 6*\dtp* are shown in 

Sec. EIH 


F. Systematic uncertainties 

We estimate systematic uncertainties in the mea¬ 
surement of the differential cross section; these are 
summarized in Table ED 


1. Uncertainties for the efficiency evaluation 

The detection efficiency is evaluated using the 
signal-MC events. However, the simulation has some 
errors or ambiguities in the reproduction of detector 
performance. This translates into uncertainties in 
the efficiency evaluation. 

Tracking has a 1% uncertainty, which is estimated 
from a study of the decays D*± —> D —► 
K° s (-> 7r + 7r )7r + 7r (0.35% per track) as well as an 

uncertainty of the radiation by an electron within 
the CDC volume. 

The electron identification efficiency in this mea¬ 
surement is very high, around 98%, and a 1% sys¬ 
tematic uncertainty is assigned for this term. Detec¬ 
tion of the 77 pair for reconstructing 7rl has a 3% 
uncertainty. The mass resolution of 7rl is well repro¬ 
duced in the MC simulation but the selection with 
the invariant mass introduces an additional 3% un¬ 
certainty, according to a comparison of the mass dis¬ 
tributions between the MC and data samples. The 
7 t 2 reconstruction has a 3% uncertainty also. The 
uncertainty of the n° reconstruction efficiency is esti¬ 
mated by a comparison of the yields of D° to 
and K^^tt 0 decays. We combine the two uncer¬ 
tainties for 7r 1 in quadrature (then 4.2% for 7rl) and 
add the uncertainty for 7t2 to it linearly, because the 
uncertainties for the two pions are fully correlated: 
the overall uncertainty is 7.2% for the two pions. 

A kinematical condition using E rat io applied in the 
selection gives an uncertainty of 2%. In addition, 
ambiguities in the detector edge locations and other 
geometrical-definition effects cause an uncertainty of 
2 %. 


The uncertainty from the trigger efficiency is es¬ 
timated by changing the energy thresholds of the 
Bhabha-veto in the HiE trigger (1%) and the clus¬ 
ter energy of the Clst4 trigger (2%), as performed in 
our previous analyses mu m The estimation 
methodology is the same but the effects are different 
among the different processes. The total uncertainty 
for the trigger efficiency is 3%. 

Background tracks and photons overlapping with 
the signal events may reduce the efficiency; this effect 
is accounted in MC by embedding the non-triggered 
event pattern (random trigger) in the signal. We 
estimate this effect by investigating the background 
conditions during different beam conditions or run 
periods. The effect on the efficiency is estimated to 
be 2%. The largest effect for the present analysis 
is estimated to come from background photons that 
form an extra n 0 with another true or background 
photon. 

We take into account an uncertainty for the 
efficiency-correction factor arising from the ip* de¬ 
pendence. Half of the difference of the factors in the 
two neighboring Q 2 regions is assigned to the uncer¬ 
tainty for each (W ,| cos0*|) region for the evaluation. 
We choose the larger for the 8 GeV 2 < Q 2 < 12 GeV 2 
regions, where two neighboring Q 2 regions exist. The 
systematic error ranges from 1% to 16%, according 
to strength of the Q 2 dependence of the factor. 


2. Uncertainties from the background subtraction 

We take into account the backgrounds from the 
single 7r° production process and the virtual Comp¬ 
ton process. We do not perform any qualitative sim¬ 
ulation for the latter process and use only an es¬ 
timate based on the features found in the experi¬ 
mental data. We assign the subtraction size applied 
for these two background sources as the systematic 
uncertainty for |cos0*| > 0.9. In addition, we as¬ 
sign another 1% error, added linearly, from the same 
sources for all the angular bins. 

We assign half of the subtraction size as the un¬ 
certainty of the background subtractions for the 77r° 
production and non-exclusive processes; they are es¬ 
timated using the signal and background MC and/or 
the experimental sideband events. 

We add the systematic uncertainty from the latter 
backgrounds in quadrature, but the others linearly, 
because the first three have a common feature that 
photon(s) and pion(s) are from backgrounds. The 
last one is different: some particles escape detection. 

The total uncertainty in the background subtrac¬ 
tion ranges from 1% to 23%, depending on the ( W , 
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| cos 0*|, Q 2 ) bin. The error is relatively large for 
the forward angles in the lowest W or the higher Q 2 
regions. 


3. Uncertainties from other sources 

The unfolding procedure has an uncertainty of 2%. 
The radiative correction has an uncertainty of 3%. 
The evaluation of the luminosity function gives an 
uncertainty of 4%, including a model uncertainty 
for the form factor of the untagged side (2%). This 
model uncertainty is based on the difference of the 
product of the luminosity function and efficiency be¬ 
tween the two models: one with ~ 1/(1 + Q 2 /m 2 ) 
and the other constant in Q\ (effectively the same 
as the model with a dependence of ~ l/(Qi + Q|) d 
because Q\ <C Q 2 , where d is a power parameter 
representing the high-Q 2 behavior). The integrated 
luminosity has a measurement uncertainty of 1.4%. 

The systematic uncertainties are added in quadra¬ 
ture unless noted above. The total systematic un¬ 
certainty is between 11% and 26%, depending on the 
{W, | cos d*|, Q 2 ) bin. 

In the following analysis, we treat a 4% systematic 
uncertainty — from the kinematical cut, geometrical 
acceptance and partially from the trigger efficiency 
and unfolding for Q 2 — a bin-by-bin error that dis¬ 
torts the W dependence. We also take into account 
the systematic uncertainty in the efficiency correc¬ 
tion from ip* dependence, which distorts the | cos0*| 
dependence. The remaining part is treated as an 
uncertainty of the overall normalization. 


VII. MEASUREMENT OF TRANSITION 
FORM FACTORS 


TABLE VI: Sources of systematic uncertainties. The val¬ 
ues indicated in a range show the range sizes in different 
bins. 


Source 

Uncertainty (%) 

Tracking 

1 

Electron-ID 

1 

Pion-pair detection (for two 7 r°’s) 

7.2 

Kinematical selection 

2 

Geometrical acceptance 

2 

Trigger efficiency 

3 

Background effect for the efficiency 

2 

ip* dependence 

1 - 16 

Background subtraction 

1 - 23 

Unfolding for Q 2 

2 

Radiative correction 

3 

Luminosity function 

4 

Integrated luminosity 

1.4 

Total 

11 - 26 




***** 
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Figure [22] shows the integrated cross section as 
a function of W for the process 7*7 -A tt 0 tt° in 
nine Q 2 bins. Peaks corresponding to the /2(1270) 
and /o(980) are evident. In this section, we extract 
the Q 2 dependence of the TFF of the /o(980) and 
those of the helicity-O, -1, and -2 components of the 
/ 2 (1270). 

A. Partial-wave amplitudes 

Helicity amplitudes in Eq. (fTT| are functions of W, 
Q 2 , and 0* but not of ip* |14| and can be expanded 
in terms of partial waves. In this channel, only even 


FIG. 22: Integrated cross section for 7*7 —> Tv°n° in nine 
Q 2 bins in GeV 2 indicated in each panel. 


angular-momentum partial waves contribute. Fur¬ 
thermore, in the energy region W < 1.5 GeV, J > 2 
partial waves (the next having J = 4) may be ne¬ 
glected so that only S and D waves need be consid¬ 
ered. Then, Eqs. © to 0 can be written as 

to = \SY$ + D 0 Y«\ 2 + \D 2 Y 2 \ 2 +2e 0 \D 1 Y}\ 2 , 
h = 2 Cl 3? (0D;|F 2 2 | - 5*F 0 ° - D* 0 Y^)D^\) , 
t 2 = —2e 0 3? (D*\Y 2 \(SY 0 ° + D 0 Y 2 °)) , (29) 
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where S is the S-wave amplitude, Do, D\, and D 2 
denote the helicity-O, -1, and -2 components of the 
D wave, respectively m, and Yj 1 are the spherical 
harmonics: 


II 

FT 

V 47T ’ 


II 

°>? 

yA(3cos 2 r-i), 


1*2 1 = ■ 

/ 1 p) 

\ — sin0* cos#* , 

V 87r 


\Y 2 1 = ■ 

/ 15 . ■« 

T sm e ■ 

(30) 


We take the absolute values of the spherical har¬ 
monics because the helicity amplitudes do not have 
ip* dependence M■ When integrated over the az¬ 
imuthal angle, the differential cross section can then 
be expressed as 


da( 7*7 —> 7r°7T 0 ) 


47rd| cos 9 
+2e 0 \D 1 Yi\" + \D 2 Y 2 


= SYo° + D 0 Y 2 ° 


0 

2 . 1 _ T r-2 12 


(31) 


The angular dependence of the cross section is 
governed by the spherical harmonics, while the W 
and Q 2 dependencies are determined by the partial 
waves. 


B. Parameterization of amplitudes 


The /o(980) and the /2(1270) have been clearly 
observed in no-tag production [28]. The sig¬ 

nal is visible in the Q 2 -integrated spectrum (Fig. 
as well as in most of the Q 2 bins in Fig. 22 Moti¬ 
vated by this, we try to extract the Q 2 dependence 
of the TFFs: Ff 2 (Q 2 ) of the /2(1270) together with 
its helicity-O, -1, and -2 components and Ffo{Q 2 ) of 
the /o(980). This is done by parameterizing S, Dq 
(and/or D{), and D 2 and by fitting the data in the 
energy region 0.7 GeV < W < 1.5 GeV. 

S and Di (i = 0,1,2) are parameterized as fol¬ 
lows (considering the general case where the in¬ 
dependent cross section is to be fitted): 

S = Af^e^+Bse 1 ^, 

A = ^r/(Q 2 )%(i2 7 o)e i ^+i?^e^-(32) 


and Bjji are non-resonant “background” amplitudes 
for the S and D waves; 4>bs, <t>BDi, 4>fo and 4>f2Di are 
the phases of background amplitudes (S and Di), the 
/o(980) and / 2 (1270) in the D waves, respectively. 
The phases are assumed to be independent of Q 2 . 
The contributions of the /2(1525) and fo(Y) (that 
are included in no-tag TT°n° production [28]) are ne¬ 
glected. The overall arbitrary phase is fixed by tak¬ 
ing 4>f2D2 = 0. When fitting the <p*-integrated cross 
section, we also set <Pbs = </bdi = 0. 

A power behavior in W is assumed for the back¬ 
ground amplitudes, which are multiplied by the 
threshold factor fj 2l+1 (l denoting the orbital angular 
momentum of the two-7r° system) with an assumed 
Q 2 dependence for all the waves: 

I3a s ( Wo/W) bs 

{Q 2 /ml + 1)<« ’ 

P 5 ap 0 ( Wp/W) bD0 
(Q 2 /m o + l) CDO 
13 5 Q 2 am ( Wp/W) bD1 
(< Q 2 /ml + l) c °i 

P 5 a D 2 ( W 0 /W ) bD2 

(Q 2 /m 2 + 1 )°m ' W 


Here, /3 = y 1 — 4m 2 0 /VF 2 is the velocity of ir° di¬ 
vided by the speed of light, m^o is the 7r° mass, 
and Wo and mo are assigned the values 1.1 GeV 
and 1.0 GeV/c 2 , respectively. Note that Adi has 
an additional factor of Q 2 to ensure that the ampli¬ 
tude becomes zero at Q 2 = 0. We set > 0 (i = 
S, Dq, Di,D 2 ) to fix the arbitrary sign of each back¬ 
ground amplitude (by absorbing the sign into their 
corresponding phases). We allow bi to have a nega¬ 
tive sign because amplitudes may still be an increas¬ 
ing function of W, but we limit \b t \ < 5; large bi 
values give a rapid W dependence, which is consid¬ 
ered unphysical. 

We use the parameterizations of the /o(980) and 
/2(1270) given in Refs. [31] [323 that are multiplied 
by their TFFs to allow a Q 2 dependence. Note that 
B(fj —¥ 7r°7r°)/H(/j —► tt+tt”) = 1/2 (because the 
fj mesons are isoscalar). For completeness, we re¬ 
produce here the parameterization of the /o(980) and 
the / 2 (1270). For the f 0 ( 980) meson, we need to take 
into account the fact that its mass is close to the KK 
threshold. The parameterization we adopt is 


B s = 
Ado = 
Adi = 
Ad 2 = 


where Aj o (g 80 ) and Af 2 ( 12 70) are the amplitudes of 
the /o(980) and f 2 ( 1270), respectively; ri(Q 2 ) is the 
fraction of the / 2 (1270)-contribution in the A wave 
with the constraints 7'o + r\ + r 2 = 1 and r, > 0; Bg 


Af 0 (980) - Ffo(Q 2 ) 



9 fo 77^/0 7T7T 1 

M% ~W 16>/37r Wo' 

(34) 
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where Ffo(Q 2 ) is the transition form factor of 
the / 0 (980), yjl + Q 2 /M 2 o is the flux factor that 
arises from definition of the luminosity function 
for a tagged two-photon cross section, fix = 

y) 1 — 4 Mx 2 /W 2 is the velocity divided by the speed 
of light for a particle A' with mass Mx in the two- 
body final states, and gf 0 xx is related to the par¬ 
tial width of the /o(980) meson via Txx(fo) = 
& 3 / 0 vx/( 1 6 ’ r %)- The factor D fo is given by the 
following expression [33]: 

D fo (W) = M 2 o -W 2 ±KH*(M /o )-H*(W) 

+3?n* (M /o ) - n* {w ), (35) 


with 


n $(W) 


Px 9f 0 xx 
167r 


i H— In 

7T 


1 — fix 
1 + fix 


(36) 


where X = f or A'. The factor fix is real in 
the region W > 2 Mx and becomes imaginary for 
W < 2 Mx- The parameter values are summarized 

in Table EII1 


TABLE VII: Parameters of the fo (980) used in the fit. 
When two errors are provided, the first is statistical, and 
the second systematic. 


Parameter 

Value 

Reference 

Mass (MeV/c 2 ) 

980 ± 20 

m 

S/o (980)71-71- (GeV) 

1.82 ± 0.03 ±g; 2 7 

m 

P 77 (keV) 

n 2Q + 0 07 

U.zy _o .06 

m 

2 / 2 
9f n KKl 9 f 0 7777 

4.21 ±0.25 ±0.21 

ESI 


We use a parameterization of the / 2 (1270) in 
Ref. [28] multiplied by the TFF, -F/ 2 (<3 2 ). The rel¬ 
ativistic Breit-Wigner resonance amplitude A R (W) 
for a spin- J resonance R of mass mx is given by 


A j r (W) 


Fr(Q 2 ) 



87t(2 J ± 1 )m R 
W 


^/T tot (W)T^(W)B(n°n°) 
m 2 R ~W 2 ^im R r tot (W) ’ 1 J 


where F R (Q 2 ) is the TFF of the resonance R and 
is the flux factor mentioned above. 
Hereinafter, we consider the case J = 2. The energy- 
dependent total width r tot (W) is given by 


Ttot(bF) = Y, T xx(W) + r ot her(W) , (38) 


TABLE VIII: Parameters of the (1270). 


Parameter 

H (1270) 

Reference 

Mass (MeV/c 2 ) 

1275.1 ± 1.2 

ESI 

Ptot (MeV) 

185.it 2 ® 

ESI 

B(nn) (%) 

84.8t?;| 

ESI 

B{KK) (%) 

4.6 ±0.4 

m 

B(w) {%) 

0.40 ±0.08 

m 

S( 7 7) (10- 6 ) 

16.4 ± 1.9 

m 

r R ((GeV/c) 1 ) 

3.62 ±0.03 

m 


where X is n, K, r], 7 , etc. For J = 2, the partial 
width r_ Y x(W / ) is parameterized as ' 2 (i 


r X x(w) 


T r B(R XX) 


qx(W 2 ) \ 5 

qx(m 2 R )J 


D 2 ( qx(W 2 )r R ) 
X D 2 (qx(m 2 R )r R ) 


(39) 


where T R is the total width at the resonance mass, 


qx(W 2 ) 

D 2 {x) 



9 + 3x 2 + x 4 ’ 


(40) 


and r R is an effective interaction radius that varies 
from 1 (GeV/c ) -1 to 7 (GeV/c ) _1 in different 
hadronic reactions 137]. For the three-body and 
other multi-body decay modes, 


Tother(bF) = T r B{R ->■ other) 


W 2 

m R 


(41) 


is used instead of Eq. ( f39j ). All parameters of the 
/2(1270) are fixed at the PDG values |34J except 
for r R , which is fixed at the value determined in 
Ref. [32]) as summarized in Table |VTH 

The normalizations of TFFs are such that 
Ff 0 ( 0) = 1.00 ± 0.11 and F /2 (0) = 1.00 ± 0.06; 
the errors reflect the uncertainties of the two-photon 
decay widths (at Q 2 = 0) of the /o(980) and 
/ 2 (1270) [313. The TFF of the / 2 (1270) for the 


helicity-i component is given by \A'i(Q 2 ).F , / 2 (Q 2 ), 
according to Eq. (32 ). 


Here, we note that there is a limitation in the 
partial-wave analysis based on Eq. ( |3~[j | (i.e., through 
the y 3 *-integrated cross section). That is, we can ex¬ 
tract information on partial waves for three out of 
the four (S, Do, Di, and D 2 ) waves only. This can 
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be understood from the fact that Eq. (31) can be 
written as a function of z = cos 9*\ as a + bz 2 + cz 4 , 
i.e., only three coefficients a,b,c are independent, 
where a , 6, c are given by combinations of the S, D 0 , 
Di, and D2 waves. This conclusion holds whether or 
not the interference term $t(S*Do) exists. 

To partially overcome this limitation, we first fit 
the y?*-dependent (but Q 2 -integrated) differential 
cross section in Sec. |V11 C| to obtain information on 
the fractions of the /2(1270) in the Do, Di, and D2 
waves. Then in Sec. |V11D| this information is used 
in the fit of the (/^-integrated cross section. 


C. Analysis of the (^‘-dependent cross section 


As described above, the ^-integrated cross sec¬ 
tion does not give information on all the partial 
waves (S, D 0 , D l5 and D 2 ). Here, we analyze the 
(^-dependent cross section to partially overcome 
this problem. Because of limited statistics, the in¬ 
dependent cross section is integrated over Q 2 ; it is 
divided into nine p* bins with a bin width of 20°, 
five equal-width |cos0*| bins, and five W bins of 
0.7 - 1.1 GeV, 1.1 - 1.2 GeV, 1.2 - 1.3 GeV, 1.3 
- 1.4 GeV, and 1.4 - 1.5 GeV. The average value of 
Q 2 (Q 2 V ) is 9.6 GeV 2 . The cross section is fitted us¬ 
ing the parameterization described above by ignoring 
the contribution of the /o(980) and Q 2 dependence. 
The values of e 0 an d ei are evaluated at Q 2 V . 

Parameters describing the assumed amplitudes are 
obtained by minimizing y 2 . To search for the true 
minimum and to identify possible multiple solutions, 
1000 sets of randomly generated initial parameters 
are employed for fits performed using M1NU1T .TDj. 
Fitted values are accepted as satisfactory solutions 
when their y 2 values are within y 2 njrl + 10 (corre¬ 
sponding to 3.2cr), where y 2 nin is the y 2 value at the 
true minimum. 

There are too many parameters to be fitted simul¬ 
taneously, particularly from the non-resonant (back¬ 
ground) amplitude that interferes with the resonant 
one in each wave. Thus we study the sensitivity of 
parameters especially in the D waves by investigating 
several sets of assumptions in the parameters. 

Two categories of fits are made: Bd 2 7^ 0 and 
Bd 2 = 0, because we consider that the interference 
between the /2(1270) and a possible non-resonant 
background is important in the D 2 wave. In each cat¬ 
egory, we try cases where both ro and r\ are non-zero 


or one of them is zero in Eq. (321. The number of so¬ 
lutions found is one or two as listed in TablellXI Here, 


solutions in which any of powers (6’s) in Eq. (331 


exceed 5 are discarded as unphysical provided that 



0 0.5 10 

Icos 0*1 


0.5 1 

Icos 0*1 


FIG. 23: |cos0*| dependence in five W bins: 0.9, 1.15, 
1.25, 1.35, and 1.45 GeV from left to right and top to 
bottom and fitted results of the ( p* nominal fit. Legend 
of lines are: solid lines show the total (black), dotted IS) 2 
(green), dashed |Do| 2 (blue), long dashed |Di| 2 (red), 
and dash-dotted |Z>21 2 (purple). 


the corresponding a parameter is not consistent with 
zero. 

Two solutions are obtained for Bd 2 7^ 0 where nei¬ 
ther ro nor rq is zero. The favored solution that has 
a smaller y 2 (referred to as a “ip* nominal fit”) has 
a small Bd 2, which naturally gives almost identical 
values with respect to the fit with Bd 2 = 0. The fit¬ 
ted values of Ff 2 (Q 2 v ) are rather similar for all the 
fits. The fitted fraction of the /2(1270) in the Do 
wave (ro) is large; the assumption of ro = 0 is dis¬ 
favored as it gives a much worse y 2 . If, in addition, 
ri is fitted, the obtained value is ri = O.I5+003 for 
the <p* nominal fit, whose value is used in the “ri 
fit” described in Sec. |VIID| In the second solution, 
Bd 2 interferes destructively with the /2(1270) giving 
a smaller value of ri = 0.11 ± 0.03. Figures [25] [27] 
show the fitted results of the ip* nominal fit that are 
projected onto the variable plotted and integrated 
over the other variables. 
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TABLE IX: Fitted parameters of the ^‘-dependent cross section. Here, = 9.6 GeV 2 . 


Parameter 


Bd 2 

7^0 



Bd 2 = 

0 


condition 

r o v o n 

n V 0 

ro = 0 

ri = 0 

r o V 0 f| r i 7^ 0 

ro 

= 0 

ri = 0 

No. 

sol. 

2 


1 

1 

1 



1 



Sol.l 

Sol.2 




Sol.l 

Sol.2 


X 2 /ndf 

236.6/208 

243.1/208 

357.4/210 

289.6/210 

241.3/211 

357.5/213 

366.0/213 

308.0/213 

Ff 2 

[Q avXxlO" 2 ) 

1.67 ± 0.15 

1.76 ± 0.11 

1 ftQ+0.05 

l.DO o.o6 

1.46 ± 0.09 

1.70 ±0.08 

1 ftQ+0-05 

I.OO-006 

1.65 ±0.06 

1.41 ±0.06 

<lD2 

(Vnb) 

0.03 ±0.03 

0-17±g;“S 

0± 0.003 

fi 90+0.02 

u - zo -0.03 


0 (fixe 

d) 


&D2 


o q+ 2.6 
**•*—8.8 

4-Oli.o 

-4.6 ±7.9 

9 0+O.6 

z ‘°-0.5 


- 



4>BD2 (°) 

89 ± 123 

114 t\l 

198 ±1 

101 ± 13 


- 



n> 


0.76 ±0.06 

0 -W-oTs 

0 (fixed) 

pj 01 +0.06 

U *®-*--0.07 

0.70 ±0.04 

0 (fi 

xed) 

0.681S:S? 

Ido 

(Vnb) 

0.13 ±0.03 

-0.13±g;g| 

0.06 ±0.02 

0.15 ±0.03 

-0.13 ±0.02 

0.06 ±0.02 

pi 90+0.02 

u ' zo -0.03 

0.16 ±0.02 

&D0 


— l-Otl.5 

-0.8l?;® 

-1.5li; 3 7 

—0.6l?i? 

-0.61?;® 

—1.511® 

1 c+0.6 
1 ’°-0.4 

-0.3 ±0.6 

4 >bdo (°) 

1861/22 

350 ± 21 

30011® 

161 ± 15 

171 ± 12 

3001?® 

271 ±9 

15111° 

0/2DO (°) 

161 ±1| 

160±ig 

- 

1441“ 

20 ± 12 



184 ±9 

ri 


0-15+oog 

0.11 ±0.03 

0 QQ+0.03 

U.UO — o 02 

0 (fixed) 

0.16 ±0.03 

n qo+0.03 

U.UO-o 02 

0.07 ±0.04 

0 (fixed) 

am 

(Vnb) 

-0.11 ±0.03 

f) 1 c+0.02 
u - lo -0.03 

-0.09 ±0.03 

0.21 ±0.01 

-0.10 ±0.02 

0.09 ±0.03 

0.15 ±0.02 

0.17 ±0.01 

bm 

l-8±l:r 

3.4t° 0 ; 9 7 

2.9ll;g 

1 Q+ 0 - 3 
-*-• J -0.2 

2.ltli 

2-911:3 

4 o+0.9 

1 q+0-3 

1 -°-0.3 

4>bdi (°) 

2 ± 5 

2251/20 

166 ±10 

209 ± 17 

140 ± 12 

345 ± 10 

356 ±8 

180 ±9 

0/2D1 (°) 

148 ± 19 

1551?? 

278 ± 29 

- 

163 ± 13 

9Q + 23 
^°-36 

69l?l 

- 

a s (Vnb) 

-0.30 ±0.03 

0.16 ±0.02 

0.29 ± 0.02 

0.16 ±0.02 

0.30 ±0.01 

0.29 ±0.02 

0.18 ±0.02 

0.32 ±0.01 

bs 


2.0 ±0.6 

—0 fi +0 ' 4 

U -°-0.5 

9 9+0.3 
Z ' Z -0.3 

-0 1+ 0 - 4 
u - i -0.4 

l-0lo.2 

9 9+O.3 
z * z -0.3 

o.elg.t 

1 4+ 0 - 2 

4>bs 

n 

262 ± 19 

70+15 
‘ °-12 

76l 6 6 

60 ±8 

95l/? 

761® 

74 ±7 

sol? 


D. Fitted results for the ^‘-integrated cross 
section 


We fit the tp*-integrated differential cross section, 
Eq. (TTJ, in the W region 0.7 GeV < W < 1.5 GeV 
for all Q 2 with Eq. ( f32j ) and Eq. ( [33] ). In the fit, the 
usual x 2 i s replaced by its corresponding quantity 
using the Poisson likelihood A defined as [3gj 


Xp = In A = 2 ^2 


Pi - rii + rii In [ — 
Pi 


(42) 


where, rii and pi are the numbers of events observed 
and predicted in the i-th bin and the sum is over 
the bins within the fitted W range. We fit the pre¬ 
dicted number of events in each bin that is related 
to the differential cross section. The latter is con¬ 
verted to the number of events by multiplying by a 
known conversion factor given from Eqs. (211 and 


(22). Here, we include zero-event bi ns in calculating 
Xpin Eq. (42). In fitting with Eq. (421, systematic 
uncertainties in the cross section are not taken into 
account. Their effects are considered separately in 
the study of systematic uncertainties in Sec. |VIIE| 
The TFFs for the /o(980) and /2(1270) and the 
fractions g(Q 2 ) are floated in each Q 2 bin, i.e., 
Ffo{Q 2 )i F f2 (Q 2 ), and r,;(Q 2 ) are obtained in each 
Q 2 bin. 


Because of the limitation mentioned above, we 
cannot determine Dq and D\ simultaneously to¬ 
gether with S and D 2 . 

Here, we prioritize the determination of D 0 and 
D 2 over D\ by setting Bp i = 0 and 


n(Q 2 ) = r 1 (Ql v ) 


Q 2 

Ql 


(43) 


where Q 2 V = 9.6 GeV 2 is the average Q 2 value for 
the Q 2 -integrated cross section and r\(Q 2 v ) = 0.15 
determined in the ip* nominal fit, and d is a free 
parameter. The effect of setting B^i = 0 is consid¬ 
ered in systematic studies. In this “ri fit,” we can 
obtain information on the helicity-1 TFF simultane¬ 
ously with those of helicity-0 and -2 in spite of the 
limitation. 

The issue with the /2(1270) is determination of its 
fractions in D 0 , D\, and D 2 in each Q 2 bin. Thus, 
in addition to the r\ fit, we perform fits assuming 
either Do V 0 with D\ = 0 (denoted as the Dq fit) 
or Pi / 0 with Do = 0 (denoted as the D\ fit). In 
each category, we put either B D2 ^ 0 or B D2 = 0. In 
the D\ fit, the evaluation of eo is necessary. It turns 
out that it has little dependence on W and cos0* 
and we use the values listed in Table [V] for Q 2 bins 
from 3.5 to 25.0 Ge V 2 . We also use the average Q 2 
value given in Table IV in each Q 2 bin in the fitting. 
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FIG. 24: p* dependence in five W bins: 0.9, 1.15, 1.25, 
1.35, and 1.45 GeV from left to right and top to bot¬ 
tom and fitted results of ip * nominal fit. Legend of lines: 
solid lines for the total (black), dotted to (red), dashed 
fi cos p* (blue), and dash-dotted f 2 COs 2 y>* (purple). 


The n fit with Bo 2 7 ^ 0 (also with Bo 2 = 0) 
gives a unique solution as summarized in Table [X] 
In the Dq fit, two solutions each are obtained for 
both Bo 2 7 ^ 0 and B 02 = 0 as also summarized 
in Table [X] In the Di fit with B 02 7 ^ 0, no so¬ 
lutions are obtained, i.e., all solutions have one or 
more powers beyond the limit. Also, in the D\ 
fit with B 02 = 0 , two solutions were found with 
Xp min /ndf = 719.7/504, which is much worse com¬ 
pared to that of the Do fit. 

As a result, it is apparent that there is a significant 
helicity-0 component of the / 2 (1270) in two-photon 
production when one of the photons is highly virtual. 
The results obtained in the ^-dependent fit that 
B 02 is non-zero and 74 is non-zero are also strongly 
supported by the p *-integrated fit. The results of 
the r 1 fit with B 02 7 ^ 0 (denoted as the “ri nominal 
fit”) give a unique solution with the minimum \p and 
are shown in Figs. [28]-[33j Figures [28]-|32] show the 
angular dependence in selected W bins in each Q 2 
bin while Fig. [33| shows the integrated cross section 
in each Q 2 bin. 

To test the contribution of the / o (980), we perform 
a fit in which the /o(980)-TFF is set to zero. This 
results in Xp /ndf = 650.2/511 to be compared to 
572.4/501 for the 74 nominal fit. The significance 
of the / o (980) contribution is 7.Id. This strongly 
supports the signature of the /o(980) and the validity 


FIG. 25: p* dependence in five |cos#*| bins: 0.1, 0.3, 
0.5, 0.7, and 0.9 from left to right and top to bottom and 
fitted results of p* nominal fit. The lines are defined in 

Kg. m 

of its TFF measurement. 


E. Estimation of systematic uncertainties of 
TFFs 

In this subsection, we estimate the systematic un¬ 
certainties for the TFFs of the /o(980) and / 2 (1270). 
They can be divided in two: the overall normaliza¬ 
tion uncertainty that affects all Q 2 bins simultane¬ 
ously and the individual uncertainties that vary in 
each Q 2 bin. The former arises from the uncertainty 
of r 77 and amounts to ± 6 % for the / 2 (1270) and ± 
11% for the / 0 (980), respectively. 

Individual uncertainties are estimated in each Q 2 
bin for the / 2 (1270)-TFF, r 0 , and / O (980)-TFF. The 
individual uncertainties considered include those re¬ 
lated to the differential cross section: normaliza¬ 
tion and distortion. The normalization uncertainty 
n in each bin of |cos0*|, W, and Q 2 varies be¬ 
tween 0.104 and 0.253; the systematic uncertainty 
from this source is estimated by multiplying the 
cross section by the factor 1 ± n. The system¬ 
atic uncertainties from distortion are estimated in 
each bin by multiplying the cross section: in W 
by the factor 1 ± 0.10(TV — 1.1 GeV) (±4% distor¬ 
tion in the full range) and in | cos 0* | by the fac- 
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TABLE X: Fitted parameters of the ip *-integrated cross section, where TFF parameters, Ff 2 {Q 2 ), ro(Q 2 ), and 
Ffo{Q 2 ), are obtained at each Q 2 bin (in GeV 2 ). 


Parameter 

ri 

fit 

Do fit 


Bdi ^ 0 

Bd 2 = 0 

Bd 2 j=- 0 

Bd2 

= 0 




Sol.l 

S01.2 

Sol.l 

Sol.2 

Xp /ndf 

572.4/501 

689.9/505 

589.2/502 591.0/502 621.2/506 622.1/506 

F f 2(0.0); (xlO~ 2 ) 



100 ±6 (def.) 



F f2 (3.5); (xlO- 2 ) 

11.5± 1.2 

i2.6ll;l 

10.9II.I 

11.6 ± 1.5 

13.611;? 

13.311;? 

F /2 (4.5);(xl0- 2 ) 

7 9+0.7 

1 -*-0.8 

7.0 ±0.4 

7.oll;l 

7.6 ±1.0 

7.5 ±0.5 

7.6 ±0.5 

F f2 (5.5); (xlO" 2 ) 

5 9+ 0 * 5 
°* y -0.6 

5.8 ±0.4 

5 5++° 
O.O-0.8 

6-Olli 

6.6 ± 0.4 

6.5 ±0.4 

Ff2(7.0); (x 10 -2 ) 

c 1 +0.3 
°* 1 -0.4 

4.8 ±0.3 

A Q+0-8 
4 *°-0.7 

5.2l?;I 

5.8 ±0.3 

5.7 ±0.3 

F f 2(9.0); (xlO- 2 ) 

4 n+0- 3 

^• u -0.4 

3.8 ±0.2 

q O+0.7 
0-O_o.6 

4 i+0- 6 
^• ± -1.2 

4.4 ± 0.2 

4.3 ±0.2 

F f2 ( 11.0); (xlCT 2 ) 

q q+0.3 
°*°-0.4 

3.0 ±0.2 

0 n+0.6 

o- u -0.5 

3.3l?;? 

3.6 ±0.3 

3.4 ±0.3 

F /2 (13.5);(xKT 2 ) 

2.5 ±0.3 

9 c+0.3 
z -°-0.2 

9 7+O.5 
z - 1 -0.4 

2.9 ±0.5 

3.4 ±0.3 

3.3 ±0.3 

F f2 ( 17.5); (xHT 2 ) 

2.3 ±0.3 

2.3 ±0.3 

2.2 ±0.4 

9 4+0.4 

Z *^-0.5 

2.9 ± 0.2 

2.8 ±0.3 

F/ 2 (25.0); (x 1CF 2 ) 

M°) 

1.8 ±0.2 
641?? 

1.4 ±0.2 

1.8 ±0.3 

11+53 

°-46 

2.0 ±0.4 

q9 + 30 
oz -89 

2.5 ±0.2 

1-615:1 

7*0(3.5); (%) 

ro(4.5);(%) 

r 0 (5.5); (%) 

ro(7.0);(%) 

ro(9.0);(%) 

r 0 (11.0);(%) 

r 0 (13.5);(%) 

r 0 (17.5);(%) 

45. It 22.1 
37.61?/? 
4i.7±it:g 
46.311?;? 
35.8+1?/ 
26.21111 
71.3 ± 19.3 
53 8+ 21,8 

00.0-21 0 

51 0 +14 ' 9 

01. u _i6.5 

2 0 +4,9 
z,u -2.0 

13.8+5? 

6 1+ 6 3 

4 . 7 + 3 / 

7 5+ 8 - 1 
10.9l?° 2 6 
26.0111:1 

37.01?/;? 

17.81/+ 7 

28.011/;/ 

24.61/+ 8 

21.41/+ 9 

91 -1+15.2 
Z1 - 1 -10.9 
1x9 c+20.0 

oz.o_ 18>1 

44.51??;? 

qq r + 16.7 
00.0 _ 14 9 

15.81/°!° 

25.61?+ 

22.5l/° 4 6 

20.3l/° 7 9 

20.7l? 3 4 6 

40.211?;? 

q 7 7+19.9 
* -14.7 

56.61??;? 

H-7151 

31.21?;/ 

31-4151 
21.81/1 
90 q+9.9 

^O.C$ — io.4 

47.911/1 

55.61?+ 

52.91??;? 

!0.3l?1 

28.81/1 

28.5+11 

18.9+/1 

22.51??;/ 

43.011?;/ 

52.61/51 

r 0 (25.0);(%) 

75.3 ±24.5 

0.1 ±4.3 

59.9I23.0 

4Q 7+21.4 
' -18.6 

62.41?°;? 

3.7 ±7.7 

<Pf 2 Do(°) 

5911? 

29611? 

359 ± 10 

348l? 8 

3001? 

304l| 

d in ?'i(Q 2 ) 

-0.2 ±0.3 

0.00 ±0.04 

- 

- 

- 

- 

F f0 ( 0.0); (xlO" 2 ) 
Ffo(3.5); (x 10 -2 ) 

18-6111 

11.81/;? 

100 ± 11 (def.) 

1 6 2 + 8 ’ 4 1 5 7 + 8 - 3 

10.91/;? 

11.51/;? 

F /0 (4.5);(xlO- 2 ) 

7.6 ±4.9 

q 0 + 2.7 
0.o_3.0 

0 7+4.0 
* -6.2 

8-31/1 

3.0 ±3.2 

3.4 ±3.2 

F/o(5.5); (x 10 -2 ) 

8.8 ±2.2 

3-011;? 

7 o+2.4 
* -0-2.6 

7.8l|? 

2 Q+ 21 

+.»_ 2 .4 

q 9 + 2.1 

°’ Z -2.4 

F /0 (7.0);(xl0- 2 ) 

7.4 ± 1.5 

3.0 ± 1.1 

7.2 ± 1.5 

7.21?;? 

3.6 ± 1.1 

3.7 ± 1.1 

Ffo(9.0); (x 10 -2 ) 

4.7 ± 1.5 

2.1 ±0.9 

5.2111 

5.H11 

2.8 ± 1.0 

2.7 ± 1.0 

F /0 (11.0);(xl0^ 2 ) 

7-4li;l 

4.0 ±0.9 

7.4 ± 1.1 

7.3 ±1.1 

4.8 ±0.9 

4.6 ±0.9 

F/o(13.5); (x 1CF 2 ) 

5.411:1 

2.3 ±0.8 

5.5 ± 1.1 

5.5111 

3.2 ±0.8 

3.0 ±0.8 

F /0 (17.5);(xl0^ 2 ) 

3.9 ± 1.0 

0.7 ±0.8 

q q + 1.0 
°*°-1.2 

3.411? 

1.8 ±0.9 

1.3 ±0.9 

Ffo(25.0); (x 10^ 2 ) 
<t>fo(°) 

O r+0.8 

°*°- 0.9 

23 ± 10 

1.0 ±0.8 

342l« 

q -1+0.9 
o - L - 1.0 

30l? o 

q q + 0.8 
O.O-0.9 

28l? 0 

1.9 ±0.8 
3561?? 

1-41/1 

3571?? 

a s (Vnb) 

9 1 + 0.5 
z * ± - 0.4 

1.8 ±0.3 

2-41/1 

2.41/1 

2.1 ± 0.3 

2.0 ±0.3 

bs 

1.1 ±0.2 

1.8 ±0.2 

0.8 ±0.2 

0 .811;? 

1.5 ±0.1 

1.5 ±0.1 

cs 

0.6 ±0.1 

0.5 ±0.1 

0.6 ±0.1 

0.5 ±0.1 

0.5 ±0.1 

0.5 ±0.1 

ciDo(Vnb) 

boo 

0.2 ±0.1 
-i-3ll:? 

0 5+ u 1 
u -°o.o 

—0 7 + 0 - 4 
u -' -0.4 

0 -211;? 
-0.3 ± 1.1 

0.211;? 

1 7+1.6 

4 - ' -4.3 

0 7+0.6 

u -' - 0.4 

9 9+1-2 
z * z —1.6 

0.4 ±0.1 
-1 6+°' 7 

Cdo 

0.0 ±1.0 

0.0 ±0.1 

0 o +0 - 4 

U.U-o 4 

0.0 ±0.8 

0.4 ±0.3 

0.0 ±0.1 

<Pdo(°) 

93l? 0 

2711/ 

57 +16 

01 -19 

95l? 7 

265l? 0 

272 ±5 

a D2 (Vnb) 

O.Slo® 

0 (fixed) 

0 5 +0 - 3 
U -°-0.2 

0.4111 

0 (fixed) 

0 (fixed) 

bD2 

0.311J 

0 (fixed) 

9 1 + 0.5 
Z “ L — 0.9 

-0.9l?;? 

0 (fixed) 

0 (fixed) 

CD 2 

0.3 ±0.3 

0 (fixed) 

0.2 ±0.2 

0.2 ±0.1 

0 (fixed) 

0 (fixed) 
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W (GeV) W (GeV) 


FIG. 26: W dependence in five |cos0*| bins: 0.1, 0.3, 
0.5, 0.7, and 0.9 from left to right and top to bottom and 
fitted results of * nominal fit. The lines are defined in 
Fig. [23] 



FIG. 27: W dependence in nine ip* bins: 10° to 170° in 
20 ° step from left to right and top to bottom and fitted 
results of <p* nominal fit. The lines are defined in Fig.|24| 


tor 1 =F [(a + 6)| cos0*| — a], where a (6) has values 
0.01 — 0.09 (0.02 — 0.16) depending on W and Q 2 
(jlj or distortion in the full range). The values 
of n, a, b are estimated in Sec. |VIF| The other in¬ 
dividual uncertainties considered are from the fitted 
W range, from the parameterization and from as¬ 


sumed constants. The fitted W ranges are changed 
to 0.65 — 1.4 GeV or to 0.75 — 1.6 GeV. The pa¬ 
rameterization of the background Bi is changed to 
W/{ 1.1 ± 0.2) GeV or to Q 2 /(1.0 ± 0.5) GeV 2 . We 


also set B\ ^ 0 to study its effect. The properties 
of the /o(980) and / 2 (1270) are changed using the 
world-average values [34.. Also tr in the /2(1270) 
parameterization is varied by its error. We have also 
taken into account the uncertainty of r\(Q 2 v ) ob¬ 
tained in the ^-dependent analysis, which is eval¬ 
uated by setting the distortion of ±4% in ip*. The 
estimated systematic uncertainty of ?t(Q 2 v ) is much 
smaller than the statistical uncertainties obtained in 
the fit. Thus, we set r\(Q 2 v ) = 0.15±0.05. A shift of 
the / 0 (980) mass ((980 ±20) MeV/c 2 ) gives rise to a 
large effect due to its proximity to the KK threshold. 
All the uncertainties are summed in quadrature in 
each Q 2 bin. For the /2(1270)-TFF, systematic un¬ 
certainties for the helicity-0, -1, and -2 components 
are calculated from the values of / 2 (1270)-TFF, ro, 
and d in rq (Eq. (431). The results are summarized 
in Table IXII 


F. Q 2 dependence of TFFs for the /2(1270) and 
/o(980) 


The obtained Q 2 dependence of TFF for the 


/ 2 ( 1270) for the rq nominal fit is shown in Fig. 34 


for the helicity-2 component. Also shown are the 
predicted Q 2 dependence by Ref. [4] and those from 
Ref. f5_ (Eqs. (]l|) and ([2])). The theoretical predic¬ 
tion of Ref. H] and of Eq. © in Ref. [5] agree well 
with data. 


Figure [35] ( [36] ) shows the helicity-0 (helicity-1) 
TFF obtained for the / 2 (1270) together with the 
prediction (4] . Note that the helicity-1 TFF is calcu¬ 
lated from i/ri(Q 2 )E/ 2 (Q 2 ), where rq(Q 2 ) is given 
by Eq. (|43| . The measured value calculated from 
r 0 at Q^~i= 0, rofO) = (3.56^;^ 2 7 f)% (28], is 
also plotted in Fig. k35] while Ref. [j] predicts zero at 
Q 2 = 0 (Table |f]. The first and second errors of ro(0) 
are statistical and systematic, respectively. The pre¬ 
diction is a factor of 1.5 — 2 larger than the measured 
values. 


Figure 37 shows the obtained Q 2 dependence of 
the TFF of the /o(980) for the rq nominal fit. 
Here, the theoretical prediction for a scalar TFF in 
Ref. [4] agrees well with the measured ones up to 
Q 2 ~ 10 GeV 2 but has steeper Q 2 dependence for 
Q 2 > 10 GeV 2 . 
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FIG. 28: Angular dependence of the cross section in the indicated W bin and results of the n fit at Q 2 = 3.5 GeV 2 
(left) and Q 2 = 4.5 GeV 2 (right). Legend of lines: solid lines for the total (black), dotted IS) 2 (green), dashed |Do| 2 
(blue), and dash-dotted |Z>21 2 (purple). 




FIG. 29: Angular dependence of the cross section in the indic ated W bin and results of the r\ fit at Q 2 = 5.5 GeV 
(left) and Q 2 = 7.0 GeV 2 (right). The lines are defined in Fig. 


28 


VIII. SUMMARY AND CONCLUSIONS 


We have measured for the first time the differential 
cross section of 7 r° pair production in single-tag two- 
photon collisions, 7*7 —> 7 r° 7 r° up to Q 2 = 30 GeV 2 
based on a data sample of 759 fb -1 collected with 
the Belle detector mm at the KEKB asymmetric- 


energy e + e collider El- The kinematical range of 
the data is 0.5 GeV < W < 2.1 GeV and | cosfT| < 
1.0 in the 7*7 center-of-mass system. 

The azimuthal angle dependence shows that the 
contribution of the helicity -0 (helicity- 1 ) component 
of the /2(1270) is large (small but non-zero). The 
differential cross section is fitted by parameterizing 
partial-wave amplitudes. The transition form factors 
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FIG. 30: Angular dependence of the cross section in the indicat ed W bin and results of the n fit at Q 2 = 9.0 GeV 2 
(left) and Q 2 = 11.0 GeV 2 (right). The lines are defined in Fig. 28 







0.5 i 

Icos 0*1 


FIG. 31: Angular dependence of the cross section in the indicat ed W bin and results of the ri fit at Q 2 = 13.5 GeV 2 
(left) and Q 2 = 17.5 GeV 2 (right). The lines are defined in Fig. |28| 


of the / 2 (1270) and / 0 (980) are measured for Q 2 up 
to 30 GeV 2 and compared with theoretical predic¬ 
tions. For the /2(1270), the helicity-0, -1, and -2 
TFFs are measured. The measured helicity-2 TFF 
of the / 2 ( 1270) agrees well with the theory predic¬ 
tion of Ref. 0] and with one of the two predictions in 
Ref. (3- The helicity-0 and -1 TFF are about a factor 
of 1.5 — 2 smaller than the prediction of Ref. [31. 


The TFF of the /o(980) is also extracted; the re¬ 
sulting Q 2 dependence agrees fairly well with the 
prediction of Ref. m for Q 2 < 10 GeV 2 but has less 
steeper Q 2 dependence for Q 2 > 10 GeV 2 . 
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TABLE XI: Resulting transition form factors for the /o(980) and /2(1270) (xlO 2 ). The first and second errors are 
statistical and systematic, respectively. 



/o(980)-TFF 


/ 2 (1270)-TFF 


Q 2 


Hcl. = 0 

Hel. = 1 

Hel. = 2 

Total 

0 

100 (def.) 

1 o q+ 0 . 6 + 21.6 

i 0 ' y —0.7—9.8 


98.2 ± 0.1^6.5 

100 (def.) 

3.45 

4.46 

5.47 

6.89 

8.92 

10.93 

13.37 

-i o '7+8-1 +5.6 

±0 - 1 - 8 .4-2.0 

7 £+4.9+4.2 

1 -°-4.9-0.6 

8.8±2.2p2 

7.4 ± i.5±S:l 

4.5 ± 1.5±l% 
7.4 ± 1.2 jlo 7 

K A+ 1 . 2 + 0 .9 
°*^-1.4-0.6 

7 7 + 2 . 0 + 2 .4 

- 2 . 1 - 0.8 
a /j+ 1 . 6 + 0.6 
^•^- 1 . 0 - 1.6 

3.8±0.8to9 

o r+0.6+0.9 
°*°-0.5-0.7 

2.4 =b 0.5±0.7 

-1 7 + 0.6 + 1 .4 
± * * -0.5-0.5 

2.1±0.4l°;| 

5.0 ±0.9to5 

q -,+0.4+0.5 
°' 1 -0.5-0.9 

2.5 ± 0.3±g;^ 
2 . 1 ± 0 . 2 ^j 

i £+ 0 . 1 + 0.2 

J--O-0.2-0.4 

1.3 ± 0.1+g'_3 
0.9 ± 0 .+q ' 2 

7.2 ±2.2l^ 

4 q+o.9+i.i 

^• y -1.5-1.6 
q q+ 0.7+0.6 
°* y -0.8-1.3 

3.2 ± 0.615 * 

9 O+0.4+0.4 
z *°- 0 .5-1.0 

9 r+0.4+0.3 
z *°- 0 .5-1.1 

0.9 ± 0.71q4 

11.5 ± 1.2±1.0 

7 9+0.7 +0.9 
‘ * z - 0 . 8 - 2.2 
c q+0.5+0.7 

O.^_o. 6 -l .2 
r 1 +0.3+0 .6 
°- 1 -0.4-0.9 

A rv+0.3+0.5 
^• u - 0 .4-1.1 
q q + 0.3+0.4 
°-°-0.4-0.8 

2.5 ±0.3l2 | 

17.23 

3.9 ± 1.0+V° 5 

1.7±0.4tgJ 

0.8 ± 0 . 1 ± 0.1 

1.3 ± 0.5±g;i 

2.3 ±0.312 2 

24.25 

q r+0.8+0.7 
0.9— 0.5 

1.6±0.3t^ 

0.6 ± 0 . 1 ± 0.1 

0.6 ±0.712:5 

1.8 ± 0.2124 



3.5 Ge\ 

jj 
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FIG. 32: Angular dependence of the cross section in 
the indicated W bin and results of the n fit at Q 2 = 


25 GeV 2 . The lines are defined in Fig. 


28 


FIG. 33: Integrated cross section in the indicated Q 2 bin 
and results of the n fit. The lines are defined in Fig. |28| 
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FIG. 34: The measured helicity-2 TFF of the /2(1270) 
as a function of Q 2 . Short (long) vertical bars indicate 
statistical (statistical and systematic combined) errors. 
The shaded area corresponds to the overall systematic 
uncertainty arising from that of r 77 . Correlations of 
uncertainties between neighboring bins exist and are in¬ 
cluded in the long vertical bars. The solid line shows the 
predicted Q 2 dependence in Table [I] by Ref. U and those 
by Ref. [5j: Eq. ([T]) (dashed line) and Eq. © (dot-dashed 
line). The theoretical curves for the helicity-2 TFF are 
normalized to 1 at Q 2 = 0. 



FIG. 35: The measured helicity-0 TFF of the /2(1270) 
as a function of Q 2 . Short (long) vertical bars indicate 
statistical (statistical and systematic combined) errors. 
The shaded area corresponds to the overall systematic 
uncertainty arising from that of T 77 . Correlations of 
uncertainties between neighboring bins exist and are in¬ 
cluded in the long vertical bars. The solid line shows the 
predicted Q 2 dependence in Table [I] by Ref. [4|. 
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FIG. 36: The measured helicity-1 TFF of the /2(1270) 
as a function of Q 2 . Short (long) vertical bars indicate 
statistical (statistical and systematic combined) errors. 
The shaded area corresponds to the overall systematic 
uncertainty arising from that of r 77 . Correlations of 
uncertainties between neighboring bins exist and are in¬ 
cluded in the long vertical bars. N ote that the helicity-1 
TFF is calculated from \Jri(Q 2 )Ff 2 (Q 2 ), where ri(Q 2 ) 
is given by Eq. ( |43| ). The solid line shows the predicted 
Q 2 dependence in Table [I] by Ref. [4J. 



FIG. 37: The measured Q 2 dependence of the TFF of 
the /o(980). Short (long) vertical bars indicate statistical 
(statistical and systematic combined) errors. The shaded 
area corresponds to the overall systematic uncertainty 
arising from that of T 77 . Correlations of uncertainties 
between neighboring bins exist and are included in the 
long vertical bars. The solid curve is the prediction for a 
scalar particle given in Table [I] by Ref. [4]. 
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